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PREFACE 


This  report  was  prepared  in  the  course  of  a  continuing  study  of 
Soviet  research  and  development  of  high-current,  high-energy  charged 
particle  beams  and  their  scientific  and  technological  applications.  The 
report  examines  Soviet  research  on  the  propagation  of  intense  relativis¬ 
tic  electron  beams  (IREBs)  through  high-pressure  air  and  gases 
(p  >  10-2  Torr),  as  reported  in  Soviet  open-source  technical  publica¬ 
tions.  A  companion  RAND  report,  Soviet  Research  on  the  Transport 
of  Intense  Relativistic  Electron  Beams  Through  Low-Pressure  Air , 
R-3309-ARPA,  August  1986,  by  the  same  author  considers  Soviet 
research  on  the  propagation  of  IREBs  through  low-pressure  air  and 
gases  (10-6  Torr  <  p  <  10”3  Torr). 

The  report  should  be  of  interest  to  specialists  in  pulsed-power  and 
high-current  charged-particle-beam  research. 

This  study  was  sponsored  by  the  Defense  Advanced  Research  Pro¬ 
jects  Agency,  Office  of  the  Secretary  of  Defense  (OSD),  and  was  carried 
out  in  the  Applied  Science  and  Technology  Program  of  the  National 
Defense  Research  Institute,  RAND's  OSD-supported  Federally  Funded 
Research  and  Development  Center. 
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SUMMARY 


The  Soviets  have  actively  pursued  research  into  the  propagation  of 
intense  relativistic  electron  beams  (IREB)  (beam  energies  above  450 
keV  and  beam  currents  in  the  tens  of  kA  range)  through  air  and  other 
gases  since  the  early  1970s.  The  present  report  analyzes  the  Soviet 
developments  of  IREB  propagation  through  background  air  at  pres¬ 
sures  from  10~2  Torr  to  atmospheric,  as  reflected  by  Soviet  open  source 
technical  literature  published  over  the  past  15  years.  A  previous  report 
by  the  same  author  considered  Soviet  developments  on  IREB  propaga¬ 
tion  through  air  at  low  pressures  (10~6  to  10~3  Torr)  [1]. 

Soviet  research  on  IREB  propagation  in  air  has  greatly  increased 
recently,  as  evidenced  by: 

•  A  large  increase  in  the  number  of  researchers  and  published 
papers  on  IREB  propagation  from  1980  to  1984. 

•  A  coincident  disappearance  from  publication  in  any  area  of  sci¬ 
ence  of  five  important  experimental  research  teams,  suggesting 
censorship  action  imposed  on  these  teams. 

The  most  important  findings  of  Soviet  research  seem  to  match  TJ.S. 
observations  regarding  propagation,  as  follows: 

1.  The  formation  of  a  plasma  channel  created  by  an  IREB  propa¬ 
gating  through  background  air  and  the  effect  of  beam  parame¬ 
ters  upon  the  plasma  channel  parameters  (and  vice  versa). 

2.  Determination  of  the  background  air  pressure  for  the  optimum 
transport  of  IREB.  Two  optimum  pressure  ranges  were  deter¬ 
mined: 

(a)  The  “ion  focused  regime” — 0.06  to  0.09  Torr 

(b)  The  “low  pressure  window” — 1  Torr 

3.  Observation  of  current  enhancement,  whereby  the  IREB-induced 
current  in  plasma  is  higher  than  the  initial  beam  current. 

4.  The  effect  of  resistive  hose  instability  upon  IREB  propagation. 

Both  the  United  States  and  the  USSR  are  concentrating  on  the 
same  effects  inherent  to  IREB  propagation  through  air.  However,  an 
important  difference  between  the  Soviet  research,  as  presented  in  the 
open  literature,  and  the  U.S.  research  is  that  the  Soviet  high  current 
experiments  are  carried  out  with  pulsed  diode  machines  in  the  1  to  1.5 
MeV  range  and  with  a  poorer  quality  beam  than  the  U.S.  experiments 
such  as  on  the  ATA  50  MeV  accelerator  and  other  beams.  There  is  a 
striking  absence  of  any  high  energy,  high  current  research  in  the 
USSR. 
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Most  Soviet  research  has  been  carried  out  at  nine  institutes,  five  of 
which  have  made  extensive  experiments  using  IREBs  with  energies  of 
0.5  to  1.5  MeV,  beam  currents  of  15  to  50  kA,  and  pulse  lengths  of  20 
to  60  ns,  with  beam  propagation  lengths  from  fractions  of  a  meter  to  3 
meters.  In  one  case,  experiments  using  beam  energies  of  2  to  20  MeV 
were  made  at  low  beam  currents. 

Soviet  theoretical  investigations  on  the  formation  of  a  plasma  chan¬ 
nel  by  the  passage  of  an  IREB  through  air  have  been  carried  out  by 
computer  simulation  made  on  the  BESM-6  computer  for  cases  of  a  1 
MeV,  30  kA  IREB,  A  new  model  used  was  more  adapted  to  higher 
pressures  and  longer  distances.  This  model  took  into  account  changes 
in  plasma  channel  parameters  due  to  IREB  and  changes  in  beam 
charge  and  current  density  profiles.  Pressures  considered  in  this  model 
were  from  1  to  100  Torr  and  the  electron  energy  inside  the  plasma 
channel  was  tens  of  eV. 

According  to  Soviet  researchers,  for  optimum  beam  transport 
through  air  it  is  imperative  to  ensure  conditions  that  allow  full  neutral¬ 
ization  of  the  IREB’s  self-fields  along  the  entire  path  of  the  beam’s 
transport.  The  propagation  is  determined  by  the  plasma  channel  con¬ 
figuration  while  the  trajectory  is  determined  by  movement  of  electrons 
that  have  lower  energy  than  the  main  component  of  the  electron  beam. 

Soviet  emphasis  on  IREB  propagation  through  fairly  high  pressure 
may  be  important  in  light  of  the  recent  increased  effort.  The  high 
pressure  range  considered  in  this  report  is  basically  applicable  to  iner¬ 
tial  confinement  fusion,  where  the  beam  must  be  transported  with 
minimum  losses  from  the  accelerator  to  the  target.  However,  in  the 
case  of  ICF,  the  gas  through  which  the  beam  propagates  usually  is  not 
air,  because  propagation  through  air  rather  than  through  other  gases 
(especially  He)  was  found  to  be  much  less  efficient.  This  pressure 
range  also  applies  to  the  passage  of  IREB  through  the  atmosphere  and 
is  pertinent  to  endoatmospheric  charged-particle  beam  weapons  and 
other  military  applications. 
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The  IREB  and  plasma  channel  parameters  being  considered  by  the 
Soviet  researchers  are  as  follows: 


IREB  Parameters 

Symbol 

Units 

Beam  energy 

Eb 

MeV,  keV 

Beam  current 

h 

kA 

Beam  current  density 

Jb 

kA/cm2 

Beam  radius 

rb 

cm 

Beam  radius  (at  input  to  drift  tube) 

rbo 

cm 

Azimuthal  magnetic  field 

Be 

G,  kG 

Time  after  pulse 

t 

ns,  ns 

Plasma  Channel  Parameters 

Background  gas  pressure 

P 

Torr 

Optimum  background  gas  pressure 

Popt 

Torr 

(for  best  IREB  propagation) 

Propagation  length 

Z 

m 

Total  length  of  drift  tube 

L 

m 

Gas  or  channel  plasma  electron  density 

ne 

cm"3 

Gas  or  channel  plasma  ion  density 

ni 

cm-3 

Characteristic  plasma  electron  energy 

tk 

eV 

Conductivity  of  plasma  channel 

a 

mho/m 

Plasma  current 

Ip 

kA 

Plasma  current  density 

JP 

kA  /cm2 

Reverse  plasma  current 

Ir 

kA 

Net  total  current 

IT 

kA 

Drift  tube  inside  radius 

Rc 

cm 

Plasma  channel  radius 

rP 

cm 

Plasma  magnetic  field  decay  time 

Td 

ns 
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I.  INTRODUCTION 


Soviet  research  on  the  propagation  of  intense  relativistic  electron 
beams  (IREB)  through  fairly  high-pressure  air  (pressure  range  0.1  to 
760  Torr)  since  the  early  1970s  included  the  study  of  the  plasma  chan¬ 
nel  created  by  the  passage  of  the  electron  beam  through  air,  the  resis¬ 
tive  hose  instability  and  its  effect  on  beam  propagation,  the  effect  of 
self-fields,  current  enhancement,  gas  expansion,  return  currents, 
inherent  beam  energy  spread,  and  other  factors. 

This  report  covers  Soviet  developments  in  IREB  propagation 
through  air  where  the  beam  is  not  focused  by  external  magnetic  fields. 
The  information  was  obtained  from  Soviet  open  publications  with 
emphasis  given  to  the  last  ten  years  of  beam  propagation  in  the  USSR. 
The  background  pressure  of  air  considered  here  varies  from  0.1  Torr  to 
atmospheric  pressure.  Soviet  work  on  IREB  transport  through  air  at 
low  pressures  (from  10-6  to  10“2  Torr)  was  considered  in  a  previous 
report  [1].  Earlier  work  on  Soviet  developments  of  IREB  propagation 
through  air  can  be  found  in  an  earlier  RAND  publication  by  Kassel 
and  Hendricks  [2], 

Most  Soviet  beam  propagation  research  has  been  carried  out  at  nine 
institutes  as  shown  in  Table  1.  Researchers  and  their  leaders  are  listed 
under  their  institutes  in  the  appendix.  The  chronology  of  the  key 
Soviet  research  on  IREB  propagation  through  air  is  presented  in  Fig.  1 
for  the  nine  institutes  and  their  leading  researchers. 

Figure  1  also  indicates  the  extent,  initiation,  and  termination  of  the 
specific  research  within  the  individual  institutes.  Note  the  short  time 
span  of  the  experimental  research  carried  out  by  the  Yu.  P.  Usov  team 
at  IYaF-TPI  where  many  researchers  apparently  made  a  concentrated 
effort,  and  the  similar  short  and  limited  effort  made  at  NIIEFA  and 
NIITP  but  with  no  recent  publications.  Also  of  interest  is  the  large 
effort  carried  out  by  two  teams  at  KhFTI,  with  a  complete  cutoff  in 
publications  in  1981  and  1982,  respectively.  With  the  termination  of 
publication  on  these  five  important  experimental  investigations  comes 
the  onset  of  five  research  efforts  that  are  currently  represented  by 
recent  publications  in  IYaF-TPI  (Team  1:  V.  P.  Grigor’yev),  IAE  (L.  I. 
Rudakov),  FIAN-IOF  (A.  A.  Rukhadze),  VEI  (A.  V.  Zharinov,  M.  A. 
Vlasov),  and  MFTI  (K.  V.  Khodatayev). 

Seven  of  the  abovementioned  institutes  carried  out  experimental 
investigations  on  IREB  propagation  through  air.  Table  2  shows  the 
accelerators  used  in  this  research,  together  with  the  beam  and  propaga¬ 
tion  parameters. 
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The  experimental  research  carried  out  in  the  first  five  institutes 
shown  in  Table  2  utilized  IREBs  with  energies  of  0.5  to  1.5  MeV,  beam 
currents  of  15  to  50  kA,  and  pulse  lengths  of  20  to  60  ns,  with  beam 
propagation  lengths  from  fractions  of  a  meter  to  three  meters.  In  the 
case  of  the  KhFTI  research,  two  experiments  are  included  where  higher 
beam  energies  were  used  but  at  low  beam  currents. 

Theoretical  work  of  any  consequence  was  carried  out  only  at  IYaF- 
TPI,  MFTI,  and  in  a  very  limited  way  at  VEI.  Out  of  the  seven  insti¬ 
tutes  that  have  performed  experiments,  only  IYaF-TPI  has  carried  out 
computer  simulation  work  of  any  magnitude.  The  theoretical  work  at 
MFTI  is  somewhat  of  a  mystery,  the  papers  not  being  clear  on  applica¬ 
tions  and  analysis.  However,  the  researchers  at  MFTI  under 
Khodatayev’s  leadership  are  of  interest,  and  their  citations  refer  to 
many  U.S.  researchers  working  on  IREB  propagation. 

Figure  2  (a  and  b)  plots  the  total  number  of  Soviet  open  source 
papers  published  specifically  on  research  of  IREB  propagation  through 
and  injection  into  air  for  the  past  15  years,  and  the  number  of  authors 
of  these  papers.  The  research  papers  on  the  use  of  IREBs  for  optical 
pumping,  beam  welding,  and  other  electron  beam  applications  that  are 
not  obvious  for  beam  transport  are  not  included.  The  number  of 
papers  as  a  function  of  years  gradually  increases  with  distinct  peaks 
during  1974,  1978,  1981,  and  1984.  The  great  jump  of  activity  in  1984 
demonstrates  an  increase  of  more  than  a  factor  of  three  of  papers  pub¬ 
lished  in  1984  over  the  number  published  in  1983. 

The  plot  of  number  of  authors  publishing  the  research  papers  is  a 
demonstration  of  relative  effort.  Each  author  was  considered  to  be 
active  (and  was  therefore  counted)  a  year  before  the  date  of  his  first 
published  paper  until  one  year  after  his  last  published  paper.  Figure  2b 
demonstrates  a  sharp  increase  in  the  number  of  authors  involved  in  the 
IREB  transport  research  from  1980  to  1984. 

This  abrupt  jump  in  the  number  of  publications  and  authors  partici¬ 
pating  in  the  research  in  1984  may  signify  an  increase  in  the  recent 
overall  Soviet  effort  on  IREB  transport  through  air,  but  it  may  also  be 
the  result  of  a  Soviet  declassification  of  earlier  research,  which  is  only 
now  able  to  be  published.  The  peak  may  therefore  represent  the 
“dumping”  of  papers. 

The  existence  of  an  active  Soviet  CPBW  program  cannot  be  directly 
inferred  from  the  available  data  as  presented  in  the  Soviet  open  source 
literature  because  of  the  low  energies  (of  1.5  MeV  and  below)  and  beam 
propagation  distances  (of  3  meters  or  less)  used  in  the  experiments. 
The  Soviet  work,  however,  does  demonstrate  a  broad  research  effort, 
involving  some  key  elements  of  IREB  propagation  through  the  air. 
The  sharp  increase  in  published  Soviet  papers  in  recent  years  implies  a 
significant  increase  in  the  research  effort  in  this  area. 


€ 


Years 

Fig,  2 — Number  of  Soviet  open  source  papers  (lower  curve)  published 
on  IREB  propagation  through  air,  and  the  number  of  authors  of 
these  papers  (upper  curve)  as  a  function  of  years 
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THEORY  AND  COMPUTER  SIMULATION 

Propagation  research  of  the  Nuclear  Physics  Institute  attached  to 
the  Tomsk  Polytechnical  Institute  (IYaF-TPI)  has  been  performed  by 
two  separate  research  teams.  The  first,  headed  by  V.  P.  Grigor'yev, 
carried  out  theoretical  calculations  and  computer  simulation  studies  of 
IREB  propagation  from  1977  to  the  present.  The  second  team,  headed 
by  Yu.  P.  Usov,  performed  several  experiments  with  IREB  propagation, 
described  in  the  next  subsection. 

In  1979-1980  the  theoretical  team  carried  out  computer  simulation 
studies  using  the  BESM-6  computer  to  investigate  plasma  channels 
formed  by  the  injection  of  IREBs  into  air  and  neutral  gas  at  pressures 
in  the  range  of  1  to  100  Torr.  The  team  considered  plasma  channel 
inhomogeneities  along  the  beam  axis,  the  effect  of  the  Bz  field  gradient 
along  the  axis  and  its  effect  upon  beam  radius,  and  beam  current  neu¬ 
tralization  in  the  plasma  channel  for  beam  pulse  lengths  of  10  to  80 
nsec.  For  these  calculations  the  researchers  took  the  values  of  the 
beam  parameters  from  early  transport  experiments  using  the  “Tonus” 
accelerator  with  a  beam  energy  of  1  MeV  and  beam  currents  of  20  to 
30  kA  [3]. 

The  computer  simulation  studies  as  carried  out  in  1979-1980  and 
described  in  Refs.  4  and  5  did  not  take  into  account  the  changes  in 
charges  and  electron  beam  current  distribution  within  the  beam  during 
IREB  transport.  The  data  and  its  discussion  as  presented  in  Figs.  3 
through  12  are  based  on  the  simpler  models,  which  are  applicable  only 
at  low  background  air  pressures  and  short  beam  transport  distances 
and  have  been  investigated  in  the  pressure  range  from  0.2  to  40  Torr. 
Computer  simulation  investigations  carried  out  in  1983  as  described  in 
Ref.  6  presented  a  more  sophisticated  model  taking  into  account  the 
beam's  charge  and  current  distribution.  These  more  recent  data  are 
presented  in  Figs.  13  through  20  and  were  determined  in  the  back¬ 
ground  air  pressure  range  from  1  to  100  Torr. 

The  parameters  involved  in  the  above  computer  simulations  were: 
the  beam  current,  J&;  total  system  current  It  [It  -  h  +  Ip>  where  Ip 
is  the  plasma  current];  beam  current  density,  J the  azimuthal  mag¬ 
netic  field  beam  radius,  r^;  electron  plasma  density,  ne\  characteris¬ 
tic  energy  of  plasma  electrons  and  the  conductivity  of  the  plasma 
channel,  <r,  as  functions  of  background  gas  pressure,  p;  propagation 
length,  z\  time  after  the  pulse,  t;  beam  radius,  r*,;  and  drift  tube  inside 
radius,  Rc.  (See  the  Glossary  for  symbol  definitions.) 
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The  variation  of  Be  and  It  along  the  beam  axis  for  background  air 
pressures  of  1  and  10  Torr  are  shown  in  Fig.  3,  and  the  axial  variation 
of  Be  as  a  function  of  p  is  shown  in  Fig.  4.  The  radial  Bd  field  distribu¬ 
tion  for  a  background  pressure  of  1  Torr  is  shown  in  Fig.  5  [4]. 


Fig.  3 — Axial  profile  of  the  azimuthal  magnetic  field  B0(solid  lines)  and 
the  total  current  IT  (dashed  lines)  for  p  =  1  and  10  Torr  for  times  t  = 
50  ns  -  curve  1,  30  ns  -  curve  2,  and  10  ns  -  curve  3  [4] 


Fig.  4— Axial  profile  of  the  azimuthal  magnetic  field  Bd  in  a  1  meter 
long  drift  tube,  t  =  30  ns  [4] 


k  A 
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1  -  .08  m 

2-  1  m 

3-  1.92  m 


Fig.  5 — Radial  B0  field  distribution  for  p  =  1  Torr,  t  *  20  ns  as  a 
function  of  various  axial  distances  [4] 


Soviet  theory  points  out  that  the  decrease  in  B$  indicates  a  notice¬ 
ably  defocusing  condition  of  the  beam  at  t  >  30  ns;  but  at  short  beam 
transport  distances,  the  basic  structure  of  the  electron  beam  does  not 
change  even  at  the  high  air  pressure  of  40  Torr. 

As  seen  in  Fig.  5,  the  magnetic  field  intensity  is  located  in  a  narrow 
band  at  the  edge  of  the  beam  with  the  area  in  the  middle  of  the  beam 
practically  magnetically  neutralized.  This  Be  field  radial  distribution  is 
the  same  for  other  background  pressures,  except  that  the  degree  of 
magnetic  neutralization  in  the  center  of  the  beam  decreases  with  an 
increase  in  pressure  (as  seen  in  Fig.  4). 

The  longitudinal  profiles  of  the  plasma  density  in  the  plasma  chan¬ 
nel  at  different  pressures  are  shown  in  Fig.  6.  The  azimuthal  magnetic 


Fig.  6 — Plasma  density  longitudinal  profile  at  p  =  1  Torr  (dashed 
line)  and  10  Torr  (solid  line)  at  times  of:  t  =  50  (1),  30  (2), 
and  10  ns  (3)  [4] 
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field  Be  decreases  along  the  beam  path  for  IREB  transport  over  dis¬ 
tances  greater  than  1  meter  in  the  pressure  range  >  10  Torr.  Such  a 
Be  gradient  demonstrates  a  significant  erosion  of  the  beam  front.  This 
effect,  however,  can  find  an  application  in  shortening  the  beam  pulse. 
The  simulation  also  showed  that  the  radius  of  the  plasma  channel  can 
be  considerably  larger  than  the  beam  radius  at  pressures  below  1  Torr. 
At  higher  pressures  the  radii  of  the  beam  and  plasma  channel  are  prac¬ 
tically  the  same.  The  beam  current  density  radial  distribution  as  a 
function  of  pressure  is  shown  in  Fig.  7  [5]. 


1  —  p*  1, 10  Torr 

2  —  p  =  0.2,  0.5  Torr 

Fig.  7 — Radial  distribution  of  the  beam  current  density  for  several 
background  air  pressures  [5] 

The  measured  and  calculated  values  of  the  total  beam  current  pulses 
are  shown  in  Fig.  8  for  a  z  =  30cm  [5].  Experimental  values  of  the 
total  current  are  shown  by  dashed  lines  demonstrating  good  agreement 
between  experiment  and  theory  for  these  cases. 

The  time  dependence  of  the  plasma  density  created  by  the  transport 
of  the  electron  beam  through  a  12  cm  radius  drift  tube  for  different 
pressures  is  shown  in  Fig.  9  [5].  The  plasma  density  produced  from 
the  collision  and  cumulative  ionization  in  the  area  of  the  beam  (r  =  0) 
increases  rapidly  and  reaches  a  maximum  value  of  about  3  x  10u  cm~3 
at  a  pressure  of  10  Torr  and  2  x  1014  cm"3  at  a  pressure  of  0.2  Torr. 

Outside  the  electron  beam  the  plasma  density  increases  only  by 
cumulative  ionization  processes  and  is  strongly  dependent  upon  the 
background  gas  pressure  increasing  by  more  than  two  tens  for  a 
decrease  in  pressure  from  10  to  0.2  Torr.  The  plasma  density  within 
the  beam  in  the  pressure  range  of  from  0.2  to  10  Torr  was  weakly 
dependent  upon  pressure  and  mainly  dependent  upon  beam  electron 
density. 
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2  -  Rc  -  6  cm 

3  -  Rc  -  1 2  cm 


p  -  10  Torr 

1  -  Rc  -  4.6  cm 

2  -  Rc  =  6  cm 

3  -  Rc  -  12  cm 

lb  -  beam  current 
dashed  line  -  experiments 
with  Rc  -  4.6  cm 


Fig.  g — Total  beam  current  pulse  as  a  function  of  background  pressure 
and  drift  tube  radius  Rc  [5] 


Fig.  9— Time  dependence  of  the  plasma  density  (ne)  as  a  function  of 
background  pressure  for  r  =  0  and  r  =  5  cm  [5] 
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The  radial  distribution  of  the  plasma  density  as  a  function  of  several 
pressures  and  several  beam  drift  tube  radii  is  shown  in  Fig.  10.  In  the 
0.2  to  10  Torr  pressure  range,  the  plasma  density  follows  the  radial 


Rc  =  4.6  cm,  1  -  p  =  0.2  Torr,  2  -  0.5  Torr,  3  -  1  Torr 
Rc  =  6  cm,  4  -  p  =  0.2  Torr,  5  -  0.5  Torr,  6-1  Torr 
Rc  =  12  cm,  7  -  p  «  0.2  Torr,  8  -  0.5  Torr,  9  -  1  Torr 

Fig.  10 — Radial  plasma  density  distribution  (ne)  as  a  function  of 
background  pressure  and  drift  tube  size  for  a  time  of  37  ns 
after  initiation  of  the  beam  pulse  [5] 


profile  of  the  beam  density,  depends  weakly  upon  pressure,  and  falls  off 
rapidly  beyond  the  beam  with  increased  background  pressure.  The 
radial  plasma  density  distribution  was  flat  at  the  beam  axis  and  steep 
at  the  beam  edge.  The  dependence  of  the  plasma  channel  radius  upon 
pressure  for  sever aljirift  tube  radii  is  shown  in  Fig.  11. 

The  characteristic  plasma  electron  energy  (c^)  averaged  over  the 
beam  cross  section,  as  a  function  of  time  for  various  background  pres¬ 
sures,  is  represented  in  Fig.  12a.  The  radial  distribution  of  the  plasma 
electron  energy  is  shown  in  Fig.  12b  for  the  various  background  pres¬ 
sures.  Figure  12a  shows  that  ek  has  a  maximum  at  the  start  of  the 
beam  current  pulse  when  the  fields  produced  at  the  beam  front  are 
large  and  then  it  decreases  with  time.  In  the  radial  distribution  of  e& 
(Fig.  12b)  the  value  of  in  the  area  of  the  electron  beam  remains 
fairly  constant  and  falls  off  slowly  with  increasing  radial  distance.  The 
value  of  €k  decreases  with  increasing  background  pressure,  which  is  due 
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Pt  Ton 

Fig.  11— Plasma  channel  radius  r  as  a  function  of  background  pressure 
for  several  drift  tube  radii  [5] 


Fig.  12 — Characteristic  average  plasma  electron  energy  (€k)  as  a 
function  of  time  (a)  and  its  radial  distribution  (b)  at  time 
t  =  37  ns  and  drift  tube  radius  of  12  cm  for  various 
background  pressures  [5] 
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to  the  increase  of  the  collision  frequency  vem  of  the  plasma  electrons 
with  the  gas  molecules  as  the  pressure  increases.  The  plasma  electrons 
lose  their  energy  with  each  collision  in  time  l/vem  and  do  not  have 
enough  time  to  obtain  a  high  energy. 

The  beam  was  found  to  induce  a  current  flowing  opposite  to  the 
direction  of  the  electron  beam  in  the  plasma.  Most  of  this  current 
flows  within  the  electron  beam;  and  its  current  density,  Jp>  was  found 
to  decrease  rapidly  with  increasing  radius. 

The  IYaF-TPI  researchers  state  that  they  are  quite  satisfied  with 
the  mathematical  model  as  presented  in  Refs.  5  and  7  that  demon¬ 
strated  good  agreement  between  experiment  and  theory  in  the  pressure 
range  of  0.2  to  10  Torr.  The  investigators  determined  an  optimum 
beam  neutralization  at  the  pressure  of  about  1  Torr  and  found  that  it 
is  very  dependent  upon  the  inner  radius  of  the  metal  tube.  They  found 
that  this  optimum  beam  neutralization  increases  with  an  increase  in 
drift  tube  radius. 

The  theoretical  models  of  Refs.  4  and  5  allowed  the  IYaF-TPI 
researchers  to  narrow  the  gap  between  theory  and  experiment.  How¬ 
ever,  in  these  models  the  beam  characteristics  were  predetermined  and 
the  change  in  charge  and  current  distributions  within  the  beam  were 
not  taken  into  account  in  the  process  of  beam  transport.  These 
assumptions,  as  was  shown  in  the  experiments  of  Refs.  3  and  7,  are 
correct  only  in  the  case  of  low  pressures  and  short  distances,  because  at 
higher  background  gas  pressures  and  larger  beam  transport  distances, 
the  beam  current  and  charge  density  profiles  change,  which  leads,  in 
turn,  to  a  change  in  the  plasma  channel  parameters.  Another  model 
for  use  at  higher  pressure  and  longer  transport  distances  was  proposed 
in  1983  in  Ref.  6,  taking  into  account  the  change  in  the  parameters  of 
the  plasma  channel  formed  by  the  IREB  and  the  changes  in  beam 
charge  and  current  density  profiles. 

A  numerical  simulation,  made  in  Ref.  6,  assumed  an  IREB  through 
air  at  pressures  from  1  to  100  Torr,  with  the  electron  energy  in  tens  of 
eV  inside  the  plasma  channel,  and  approximating  the  recombination 
coefficient  by  the  Mehr  and  Biondi  formula  [8].  The  plasma  density 
ne  as  a  function  of  transport  distance  and  time  is  shown  in  Figs.  13 
and  14,  respectively  [6]. 

The  plasma  density  was  seen  to  reach  4.5  x  1015  cm-3  in  30  ns  at 
pressures  50  to  70  Torr  and  about  2.5  x  1016  cm“3  at  pressures  from  1 
to  20  Torr.  Beyond  30  ns  the  plasma  density  decreased  because  the 
recombination  rate  (mostly  dissociative)  was  greater  than  the  ioniza¬ 
tion  rate.  The  plasma  channel  radius  at  pressures  greater  than  1  Torr 
coincided  with  the  beam  radius,  and  the  radial  profile  ne(r,z,t)  was 
very  close  to  rib(r9z,t).  The  plasma  channel  was  inhomogeneous  along 
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ntlr=0)t1CT,s,  cm'5 


Dashed  lines  represent  a  pressure  of  1  Torr 

1  -  t  =  10  ns  2  -  t  *  40  ns 

Solid  lines  represent  a  pressure  of  20  Torr 

1  -  t  =  6  ns  2  -  t  =  10  ns 

3  ~  t  =  17  ns  4  -  t  =  40  ns 

where  rb0  is  the  beam  radius  at  the  entrance  to  the  drift 

tube  and  Rc  is  the  drift  tube  radius 


Fig.  13 — The  plasma  electron  density  (ne)  distribution  along  the 
drift  tube  with  rbo  =  1  cm,  Rc  -  2  cm  [6] 

the  length  of  the  drift  tube,  especially  at  the  higher  pressures  and  for 
large  initial  beam  densities. 

The  dependence  of  the  density  of  the  plasma  in  the  channel  upon 
gas  pressure  and  initial  beam  current  density  was  found  to  be  complex, 
because  of  all  these  factors  affecting  beam  density  and  ionization  rate. 
At  short  distances  from  the  point  of  beam  injection  the  plasma  electron 
density  ne  was  greater  at  higher  gas  pressures.  The  ne  dependence  was 
different  at  the  output  of  the  drift  tube,  however.  For  t  <  10  to  12  ns, 
ne  was  lower  at  p  «  1  Torr  than  at  p  =  10  to  100  Torr.  However,  for 
t  >  12  ns,  ne  increased  more  rapidly  at  p  -  1  Torr. 

The  characteristic  energy  of  the  plasma  electrons  (e*)  was  found  to 
depend  only  weakly  upon  the  gas  pressure.  The  value  of  e*  was  found 
to  vary  from  1  to  14  eV  over  this  pressure  range  and  along  the  whole 
length  of  the  channel.  The  higher  values  of  €*  (5  to  14  eV)  were  found 
for  p  =  1  Torr  and  the  lower  values  (0.5  to  2  eV)  for  p  =  100  Torr. 
The  value  of  was  found  to  be  very  weakly  dependent  upon  the 
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n,(r=0)*f0'r5t  cm* 


Solid  lines  represent  Z  =  8  cm,  with  pressures  of: 

1  -  1  Torr  2-10Torr  3 -20  Torr 
4-50  Torr  5-70  Torr 

with  rbo  =  1  cm,  Rc  =  6  cm 
6-20  Torr,  rbo  =  3  cm,  Rc  =  6  cm 
Dashed  lines  represent  Z  =  92  cm,  with  pressures  of: 
1-1  Torr  2-10  Torr  3  -  20  to  70  Torr 
with  rb0  =  1  cm,  Rc  =  2  cm 
4-20  Torr,  rb0  =  3  cm,  Rc  =  6  cm 


Fig.  14 — Dependence  of  ne  upon  time  [6] 


radius.  This  can  also  be  seen  in  Fig.  12  with  earlier  calculations  of  e* 
as  a  function  of  r  and  t  [5]. 

The  plasma  channel  conductivity  (<x)  was  found  to  have  a  strong 
dependence  upon  gas  pressure  and  at  pressures  greater  than  1  Torr  was 
determined  totally  by  plasma  electron  collisions  with  neutrals.  This  a 
depends  not  only  upon  ek  but  also  upon  the  ratio  of  the  densities  of  the 
plasma  and  neutral  molecules.  The  conductivity  was  thus  found  to 
decrease  rapidly  with  an  increase  in  gas  pressure,  and  its  radial  distri¬ 
bution  is  similar  to  that  of  the  ne  distribution.  The  dependence  of 
and  <r  upon  gas  pressure  is  demonstrated  in  Fig.  15. 

The  simulation  of  the  total  current  It  -  h  +  Ip  and  the  azimuthal 
magnetic  field  Be  demonstrates  that  the  basic  plasma  current  flows  in 
the  area  of  the  axis  of  the  beam  and  decreases  as  r  approaches  r^. 
With  increase  in  gas  pressure  and  because  of  the  decrease  in  the 
plasma  channel  conductivity,  the  current  neutralization  weakens  and 
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cr,  mho/m  ek,  eV 


Solid  curves  represent  a  and  dashed  curves  represent  «k 

1  -  t  =  10  ns,  z  =  8  cm 

2  -  t  =  33  ns,  z  =  8  cm 

3  -  t  =  1 0  ns,  z  =  92  cm 

4  -  t  »  33  ns,  z  =  92  cm 


Fig.  15 — Dependence  of  the  characteristic  energy  of  the  plasma 
electrons  (ek)  and  the  conductivity  of  the  plasma  channel  (tr) 
upon  gas  pressure  (p)  for  r  =  o,  rbo  =  1  cm,  Rc  =  2  cm  [6] 


Be  grows  near  the  point  of  beam  injection.  The  radial  distribution  of 
current  neutralization  becomes  more  homogeneous.  The  Be  radial  dis¬ 
tribution  found  in  Ref.  6  is  similar  to  the  earlier  findings  made  in  Refs. 
4  and  5.  According  to  Ref.  6,  Be  increases  with  time,  and  the  sharp 
dependence  of  Be  as  a  function  of  z  is  seen  only  at  the  beginning  of  the 
drift  tube.  At  the  end  of  the  tube  the  dependence  on  t  becomes  insig¬ 
nificant,  as  demonstrated  in  Fig.  16. 

In  the  light  of  the  above  theoretical  results,  the  IYaF-TPI  research¬ 
ers  explain  electron  beam  focusing  as  follows:  A  drift  of  plasma  elec¬ 
trons  is  observed  along  the  beam  radius  as  a  result  of  the  movement  of 
the  electron  beam  in  the  Be ,  which  increases  in  time  and  falls  off  with 
distance.  Because  Be  increases  with  time,  the  electrons  obtain  a  veloc¬ 
ity  directed  toward  the  beam  axis,  and  because  of  the  Be  gradient  along 
z ,  the  velocity  is  directed  along  the  axis.  Beam  electrons  follow  trajec¬ 
tories  that  are  oscillations  with  either  decreasing  or  increasing  radial 
amplitudes,  depending  upon  which  forces  of  the  magnetic  field  predom¬ 
inate.  Thus  the  beam  either  focuses  or  defocuses. 
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BfI  kG 


Dashed  curves:  p  =  1  Torr 
1  -  t  =  7  ns 
2-t  =  17  ns 
3  -  t  =  33  ns 
Solid  curves:  p  =  50  Torr 

1  -  t  =  6  ns 

2  -  t  =  10  ns 

3  -  t  =  1 7  ns 

4  -  t  =  33  ns 


Fig.  16 — Axial  profile  of  the  azimuthal  magnetic  field  B0  for  r  =  r^  = 

1  cm,  Rc  =  2  cm  [6] 


The  beam  dynamics  calculations  are  shown  in  Figs.  17-19  [6].  In 
Fig.  17,  curve  1  represents  the  beam  pulse  without  any  current  loss  at 
the  output  with  the  following  two  conditions: 

p  =  1  Torr,  77,0  =  1  cm,  Rc  =  2  cm 

p  =  20  Torr,  7^=3  cm,  Rc  =  6  cm 

Figure  17  indicates  that  a  beam  having  a  high  current  density  of  about 
7.6  kA/cm2  at  the  input  to  the  drift  tube  can  be  transported  without 
losses  only  for  pressures  less  than  5  Torr.  At  higher  pressures,  only 
the  leading  edge  of  the  beam  pulse  can  pass,  while  the  trailing  edge  is 
cut  off.  This  shortens  the  pulse  and  decreases  its  amplitude  at  the 
drift  tube  output.  The  cutoff  effect  is  more  evident  at  higher  gas  pres-  I 
sures  and  beam  current  densities,  in  longer  drift  tubes,  and  with  1 
smaller  tube  radii.  As  an  example,  beam  transport  was  computed 
for  beam  pulse  increasing  from  It  max(z  =  0)  =  23.7  kA,  Jb0  = 

7.6  kA/cm2  to  It  max  (2=0)  =  50  kA,  Jto  =  16  kA/cm2.  Pulse  length 
was  found  to  decrease  from  33  to  26  ns  for  L  =  1  m,  Rc  =  2  cm.  The 
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It /lb,  ma* 


1  -  at  the  input  to  the  drift  tube 

2  -  7  at  the  output  of  the  drift  tube 

2  -  p  =  5  Torr 

3  -  p  «  10  Torr 

4  -  p  =  20  Torr 

5  -  p  =  30  Torr 

6  -  p  =  50  Torr 

7  -  p  a  70  to  100  Torr 

with  *b  max  =  23  7  kA-  rbo  =  1  cm,  Rc  =  2  cm 

Fig.  17 — Electron  beam  current  pulse  shape  [6] 


ratios  of  the  beam  current  values  at  the  input  and  output  of  the  drift 
tube  were  0.56  and  0.38,  respectively. 

The  above  theoretical  results  were  in  good  agreement  with  experi¬ 
ments  performed  in  1967  [9]  when  it  was  determined  that  current  was 
cut  off  at  the  trailing  edge  of  the  pulse  in  drift  tubes  1  meter  long  at 
pressures  above  6  Torr. 

Figure  17  indicates  that  the  current  cutoff  stabilized  at  pressures 
from  50  to  70  Torr,  and  there  was  no  further  pulse  shortening  beyond 
70  Torr.  There  also  was  no  current  cutoff  at  pressures  below  5  Torr. 
Both  effects  were  attributed  to  current  neutralization.  At  low  pres¬ 
sures,  the  gradient  of  the  axial  magnetic  field  is  flattened  out  (see  Fig. 
16);  and  therefore  there  is  no  pulse  current  cutoff  present.  The  con¬ 
ductivity  of  the  plasma  channel  is  then  high  enough  (Fig.  15)  and  the 
level  of  current  neutralization  is  also  sufficiently  high  during  the  entire 
pulse.  As  pressure  increased  beyond  5  Torr,  plasma  channel  conduc- 
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rb»  CM 


(a)  p  =  1  Torr;  1  -  t  =  1 0  ns,  2  - 1  =  33  ns, 

3  -  t  =  65  ns 

(b)  solid  curves  -  p  =  50  Torr:  1  -  t  =  10  ns, 

2  -  t  =  17  ns,  3  -  t  =  33  ns 
dashed  curves  -  p  =  70  to  100  Torr: 

1  -  t  =  10  ns,  2  -  t  =  33  ns 

(c)  p  =  20  Torr:  1  -  t  =  10  ns,  2  - 1  =  17  ns, 

3  -  t  =  33  ns,  4  -  t  =  65  ns 


Fig.  18 — Dependence  of  the  beam  radius  rb  upon  the  axial  length  z  [6] 


tivity  decreased  (Fig.  15);  and  the  level  of  current  neutralization  was  no 
longer  sufficient  to  flatten  out  the  Bq  field  gradient  and  to  ensure  the 
beam  transport  without  loss.  Above  50  Torr,  the  conductivity  of  the 
channel  decreased  to  such  an  extent  that  plasma  current  had  practi¬ 
cally  no  effect  upon  the  axial  magnetic  field  distribution.  Be  was  then 
determined  only  by  the  inhomogeneity  of  the  injected  beam  at  the 
beam  front  and  by  its  current  density  change  in  the  process  of  beam 
transport. 

The  change  in  beam  parameters  along  the  drift  tube  is  shown  in 
Figs.  18  and  19.  At  short  distances  from  the  beam  injection  plane 
(z  ~  8  to  16  cm),  the  beam  was  focused  by  a  sharp  increase  of  Be  in 
time.  As  the  beam  traveled  further,  it  started  spreading;  and  the  basic 
electron  losses  occurred  at  z  >  LA 2.  At  the  end  of  the  drift  tube,  at 
pressures  greater  than  10  Torr  and  T  >  18  ns,  the  electron  beam  den- 
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Rc  =  2  cm,  p  =  50  Torr_ 

1  -t  =  10ns,  2-t»14ns,  3- t  =  17  ns, 
4  -  t  =  20  ns,  5  -  t  *  26  ns,  6  -  t  =  33  ns 


Pig.  19— Relative  beam  current  losses  for  rbo  =  1  cm  [6] 


sity  on  axis  sharply  dropped,  and  the  radial  beam  profile  na(r)  was 
homogeneous  along  the  entire  cross  section  of  the  drift  tube.  Beam 
spreading  decreased  with  increasing  drift  tube  radius  Rc  (Fig.  20).  This 
can  be  attributed  to  the  increase  of  the  inductance  of  the  system  with 
increasing  Rc  and  the  increase  in  current  neutralization. 


H>,cm 


1  -  Rc  =  2  cm,  2  -  Rc  =  4  cm,  3  -  Rc  =  10  cm 
,  dashed  line  -  Rc  =  20  cm 

Fig.  20— Beam  radius  rb  as  a  function  of  distance  and  drift 
tube  radius  Rc  [6] 
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Soviet  researchers  stress  that  even  dense  electron  beams  have  little 
beam  loss  at  short  distances  from  the  beam  injection  plane  up  to  pres¬ 
sures  above  50  Torr  (Fig.  19).  A  detailed  analysis  of  the  propagation  of 
an  IREB  along  a  plasma  channel  with  a  return  current  (current  flow¬ 
ing  opposite  to  the  IREB  in  the  channel)  was  made  in  Ref.  10.  Here 
the  radial  distribution  of  the  charge  density  and  currents  of  the  elec¬ 
tron  beam,  plasma  electrons,  and  ions  have  been  determined  for  the 
beam  radius  smaller  than,  equal  to,  and  larger  than  the  radius  of  the 
plasma  channel. 

The  beam  is  assumed  to  be  in  equilibrium,  with  its  space  charge  par¬ 
tially  neutralized  by  the  background  gas  ions.  This  steady  state  condi¬ 
tion  is  possible  when  an  IREB  in  injected  into  low-density  plasma 
(where  rib  is  the  density  of  the  electron  beam  and  n;  is  the  density  of 
ions  of  the  background  gas).  However,  when  the  electron  beam  is 
injected  into  a  dense  plasma,  space  charge  neutralization  is  also  accom¬ 
panied  by  a  partial  beam  current  neutralization,  characterized  by 
plasma  electrons  flowing  opposite  to  the  direction  of  the  electron  beam. 

The  duration  of  the  return  plasma  current  flow  was  found  to  be  0.1 
microseconds  for  -  1cm,  ne  ~  1014  to  1015,  and  electron  plasma  tem¬ 
perature  of  ~  1  eV  and  is  longer  than  the  usual  beam  current  pulse 
length  [10],  As  a  result,  when  Soviet  researchers  considered  the  equi¬ 
librium  condition  of  electron  beam  transport,  they  found  it  necessary 
to  take  into  account  this  plasma  electron  current  as  well  as  the  effect 
of  the  background  gas  ions. 

The  radial  potential  distribution  as  a  function  of  the  system  configu¬ 
ration,  the  electron  beam  radius,  the  plasma  radius,  and  the  total  beam 
current  have  been  investigated  using  the  equations  developed  in  Ref. 
10.  A  specific  calculation  was  made  for  the  case  when  the  electron  and 
ion  plasma  radii  were  the  same  (r;  =  re  =  rp)  and  when  the  IREB 
radius  was  equal  to  the  plasma  radius  ( rh  =  rp). 

With  a  net  negative  charge  excess  along  the  beam  axis,  the  beam 
and  plasma  electron  density  increased  to  the  outer  beam  periphery 
while  the  ion  density  decreased  [10].  However, .  with  a  net  positive 
charge  excess  along  the  axis,  the  beam  and  plasma  electron  density 
decreased  to  the  beam  periphery,  while  ion  density  increased.  When 
the  electron  beam  current  density  on  the  axis  was  greater  than  the 
plasma  electron  current  density,  the  electron  beam  current  density 
increased  with  an  increase  in  beam  radius  while  the  plasma  current 
density  decreased.  Thus,  the  radial  distribution  of  the  total  current 
density  of  the  system  represented  a  hollow  beam  (Fig.  21a),  which  was 
in  agreement  with  the  1971  experiments  with  IREB  injection  into 
dense  plasma  [11].  The  current  transported  by  the  ions  could  be 
neglected  in  that  case. 
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(a)  (b)  (c) 


a  -  rb  =  rP>  b  “  rb  <  rP- c  “  rb  >  rP 


Fig.  21— Radial  distribution  of  the  current  density  of  the  beam,  plasma, 
and  total  current  flow  [10] 


When  the  IREB  radius  was  smaller  than  the  plasma  radius 
( rp  >  rb)y  the  current  density  of  plasma  electrons  dropped  to  zero;  with 
increasing  radius  and  changing  direction,  it  increased  monotonically  as 
it  approached  the  plasma  channel  boundary  (Fig.  21b).  The  radial  pro¬ 
file  of  the  total  current  in  this  configuration  formed  a  double  wall  tube. 

A  beam  radius  can  be  larger  than  a  plasma  radius  ( rb  >  rp)  when  the 
beam  is  transported  along  thin  metal  wires  or  along  a  thin  laser  beam. 
The  electron  beam  current  density  was  found  in  this  case  to  increase 
slightly  to  the  plasma  boundary;  beyond  it  increased  more  rapidly  to 
the  beam  boundary  as  shown  in  Fig.  21c. 

The  Soviet  researchers  state  that  IREB  propagation  in  plasma  with 
return  current  permits  effective  beam  transport  with  beam  currents 
that  are  considerably  higher  than  those  in  static  cases  of  neutralized 
electron  beams.  These  beams  have  an  application  where  it  is  necessary 
to  create  IREB  with  very  high  currents  [10]. 

IYAF-TPI  RESEARCH— EXPERIMENTS  ON  CHANNEL 
FORMATION  AND  HOSE  INSTABILITIES 

The  IYaF-TPI  experimental  research  team  headed  by  Yu.  P.  Usov 
performed  several  experiments  in  1975-1977  involving  plasma  channels 
induced  by  IREB,  hose  instability,  the  optimal  background  pressure  for 
IREB  propagation,  and  the  energy  spread  of  beam  electrons.  These 
experiments  were  carried  out  using  the  Tonus  accelerator  operating  at 
1  MeV  with  beam  currents  from  35  to  40  kA,  and  a  pulse  length  of  60 
ns.  No  publications  by  this  team  were  found  after  1977. 

The  formation  of  the  plasma  channel  was  first  observed  by  Soviet 
researchers  in  1972  when  an  IREB  was  injected  into  air  and  gases  at 
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pressures  from  1  to  760  Torr  [12].  In  1977,  IYaF-TPI  researchers 
demonstrated  experimentally  that  the  radius  of  the  plasma  channel 
created  by  an  IREB  was  greater  than  the  beam  radius  in  the  pressure 
range  of  0.1  to  1  Torr,  and  the  difference  between  the  channel  and 
beam  radii  increased  with  a  decrease  in  pressure  [13].  The  same  Tonus 
accelerator  was  used  with  a  4  cm  diameter  beam  injected  into  a  9  cm 
diameter  drift  tube  60  cm  long.  Beam  transport  was  found  to  be  most 
effective  in  the  pressure  range  of  0.09  to  0.06  Torr;  however,  long  wave 
hose  instabilities  were  observed  in  this  pressure  range  with  lateral 
beam  displacement  of  the  order  of  one  beam  diameter  (Fig.  22a). 

With  further  increase  in  background  pressure,  beam  refocusing 
increased  in  these  experiments  because  of  increased  current  neutraliza¬ 
tion.  In  the  pressure  range  for  optimum  beam  transport,  a  glowing 
cylindrical  layer  appeared  around  the  periphery  of  the  beam  (Fig.  22b). 
This  layer  was  assumed  to  be  a  plasma  channel  carrying  part  of  the 
return  current  of  the  beam.  The  channel  was  found  to  narrow  down 
with  increased  pressure;  and  when  the  beam  developed  hose  instability, 
the  cylindrical  symmetry  of  the  layer  was  not  disturbed  but  the  whole 
beam  was  simply  displaced  laterally.  Measurements  revealed  that  no 
high  energy  electrons  moving  parallel  to  the  beam  flow  were  present  in 
this  layer  region.  From  the  analysis  of  data  obtained  in  Refs.  3,  13, 
and  14,  the  IYaF-TPI  researchers  deduced  that  an  IREB  can  be  trans¬ 
ported  in  large  volumes,  in  large  diameter  insulating  channels,  and 


Fig.  22 — Photographs  of  background  gas  glow  [13] 
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through  free  space  by  forming  two  channels  for  the  return  current 
flow — one  in  the  area  of  the  electron  beam  itself  and  the  other  in  the 
thin  cylindrical  layer  around  the  beam.  The  latter  shorted  the  beam 
current  to  the  accelerator  while  maintaining  beam  self-focusing.  In 
addition,  the  plasma  layer  had  a  stabilizing  effect  upon  the  lateral 
oscillations  of  the  beam  improving  IREB  transport. 

A  quickly  developing  instability  was  also  observed  when  IREB  prop¬ 
agated  in  background  gas  at  the  pressure  of  about  100  Torr,  leading  to 
a  collapse  of  the  electron  beam  (Fig.  22c). 

Experimental  measurements  of  the  energy  spectrum  were  made  for  a 
1  MeV,  37  kA  IREB  injected  into  a  2  m  drift  tube  filled  with  air  at 
various  pressures  (Fig.  23)  [3,15,16].  Figure  23  shows  that  the  beam 
electrons  with  less  than  the  maximum  energy  are  located  in  two 
groups,  one  at  the  front  of  the  beam  pulse  and  the  other  at  the  end  of 
the  pulse.  The  amplitude  of  the  electron  current  pulse  increased  with 
increasing  electron  energy  and  peaked  at  the  maximum  electron  energy 
[3].  This  beam  energy  structure  maintained  its  form  over  the  2  m 
transport  distance  in  the  pressure  range  of  0.6  to  1  Torr.  At  20  Torr, 
the  energy  spectrum  changed  considerably  after  1  m  distance  (Fig. 
23b).  In  this  case,  both  the  energy  peak  and  the  peak  at  the  end  of  the 
pulse  are  decreased  considerably.  With  further  increase  in  background 
pressure  to  50  Torr,  both  the  energy  peak  and  the  peak  at  the  end  of 
the  pulse  have  disappeared  completely  (see  Fig.  23c).  Further  measure¬ 
ments  demonstrated  the  shortening  of  the  current  pulse  length  at  pres¬ 
sures  greater  than  20  Torr  (Fig.  24)  [3].  This  pulse  shortening  was 
attributed  to  current  cutoff  allowing  only  the  front  of  the  beam  pulse 
to  pass  through  the  dense  gas. 

It  was  also  determined  that  at  long  beam  transport  distances  the 
pulse  current  cutoff  appears  at  lower  pressures  in  the  drift  tube.  This 
current  cutoff  is  attributed  to  beam  instability  in  dense  plasma,  which, 
in  turn,  is  created  by  the  ionization  of  the  neutral  background  gas. 
The  IYaF-TPI  researchers  mentioned  in  1975  [3]  that  experiments 
were  being  carried  out  at  that  time  to  clarify  the  mechanism  and  the 
specific  instability  characteristics  that  are  responsible  for  the  above 
current  cutoff. 

The  1977  experiments  [15]  also  showed  the  existence  of  an  optimum 
pressure  of  about  1  Torr  when  the  IREB  energy  spectrum  was  least 
perturbed  during  transport.  The  maximum  energy  Em  and  its  width  at 
half  maximum  A E/Em  were  found  to  be  Em  ~  1  MeV  and 
A E/Em  ^  20  percent.  Typical  oscillograms  of  the  electron  beam 
current  pulse  as  obtained  at  the  end  of  0.8  and  2  m  drift  tubes  at  vari¬ 
ous  pressures  are  shown  in  Fig.  25. 

Beam  transport  was  found  to  be  most  efficient  when  the  beam  self- 
fields  were  fully  neutralized  along  the  whole  length  of  the  beam  path. 
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(b) 


(c) 


p  =  20  Torr 


p  =  50  Ton- 


Observed  energies: 

1— 535  to  600  keV  4—805  to  870  keV 

2— 625  to  690  keV  5—895  to  960  keV 

3— 705  to  780  keV  6—985  to  1 045  keV 


Fig.  23 — Energy  spectrum  time  scan  of  a  1  MeV  maximum  energy 
IREB  with  a  37  kA  current  pulse  at  a  1  m  transport  distance  at 
different  background  pressures  [3] 


This  occurred  at  1  Torr,  with  the  form  of  the  IREB  pulse  practically 
unchanged  (Fig.  25a  and  b,  curve  #3).  The  behavior  of  Em  as  a  func¬ 
tion  of  pressure  and  transport  length  is  shown  in  Fig.  26, 

Away  from  1  Torr,  Em  was  observed  to  be  displaced  toward  smaller 
energies  while  the  beam  energy  spread  increased.  This  behavior  of  Em 
is  not  attributed  to  energy  loss  from  gas  ionization. 

Below  popt,  the  leading  edge  electrons  were  scattered  by  space-charge 
forces.  The  nose  erosion  placed  the  Em  electrons  at  the  leading  edge  of 
the  pulse,  and  not  at  the  location  of  maximum  current  as  in  the  case  of 
p0pt  ~  1  Torr.  The  current  peak  was  displaced  toward  energies 
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(a) 


p  =  0,6  Ton 


(b) 


p  =  20  Torr 


(c) 


p  =  50  Ton- 


Fig.  24— Beam  current  pulse  oscillograph  trace  for  IREB  transport  in 
air  at  a  distance  of  1  meter  [3]  (reference  oscillation 
frequency  =  100  MHz) 


lower  than  Em,  thus  forming  the  trailing  edge  of  the  pulse  during  the 
acceleration  process.  This  was  the  reason  for  the  absence  of  the  high 
electron  energy  component  in  the  beam  spectra  and  the  Em  shift 
toward  lower  energies.  Above  popt>  the  front  part  of  the  beam  prop¬ 
agated  more  efficiently.  Above  5  Torr  (Fig.  25b),  beam  transport  ter¬ 
minated  before  the  Em  electrons  could  appear.  These  factors  lead  to 
the  same  spectra  change  as  in  the  case  for  pressures  lower  than  popt . 

The  Soviet  researchers  concluded  that  the  observed  variation  in 
beam  energy  spectra  in  beams  transported  in  neutral  gas  at  various 
pressures  is  due  to  the  different  propagation  effectiveness  of  the  vari¬ 
ous  parts  of  the  beam.  This  effect  is  maintained  by  the  transient 
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1  -  p  =  0.05  Torr 

2  -  p  =  0.1  Torr 

3  -  p  =  1  Torr 

4  -  p  =  5  Torr 


5  -  p  =  10  Torr 

6  -  p  =  15  Torr 

7  -  p  =  20  Torr 

8  -  p  =  50  Torr 


Fig.  25— IREB  current  pulses  at  the  end  of  the  drift  tube  (a)  0.8  m  and 
(b)  2  m  at  various  background  gas  pressures  [15].  The  vertical  line 
denotes  the  location  of  electrons  with  energy  Em 


IREB  energy  structure  determined  by  the  acceleration  process.  To 
maintain  an  unchanged  energy  character  of  the  beam  under  transport 
it  is  necessary  to  insure  conditions  that  allow  full  neutralization  of  the 
beam’s  self- fields  along  the  whole  path  of  transport  [15]. 

Control  of  the  length  of  the  IREB  current  pulse  was  investigated  in 
1977  using  the  principle  of  beam  current  cutoff  during  its  transport  in 
jueutral  gas  at  pressures  higher  than  popt  [17],  These  tests  were  also 
carried  out  with  the  Tonus  accelerator  at  1  MeV,  40  kA,  and  a  drift 
tube  2  meters  long.  The  current  pulse  length  was  controlled  by  varying 
the  pressure  within  the  drift  tube.  The  dependence  of  the  pulse  length 
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Em,keV 


1  -  L  =  0.3  m 

2  -  L  =  0.8  m 

3  -  L  =  1 .3  m 

4  - L  =  2  m 

Fig.  26 — Value  of  the  energy  at  the  maximum  current  peak  Em  as  a 
function  of  pressure  p  and  transport  length  L  [15] 


on  background  pressure  and  transport  length  is  shown  in  Fig.  27. 
Pulses  shorter  than  10  ns  can  easily  be  obtained  by  this  method. 

Other  experiments  at  IYaF-TPI  made  in  1977  on  plasma  channels 
induced  by  an  IREB  propagating  through  gas  and  air  indicated  that  the 
beam  tracked  the  plasma  channel  configuration.  These  experiments 
were  made  with  the  Tonus  accelerator  with  a  1  MeV,  36  kA,  60  ns 
beam  injected  into  a  drift  space  filled  with  air  [14].  The  most  effective 
IREB  transport  was  achieved  at  a  background  air  pressure  of  0.09  to 
0.6  Torr.  At  the  pressure  range  of  0.14  to  0.17  Torr  long  wavelength 
“hose”  instabilities  were  observed  with  a  large  lateral  displacement  of 
the  beam.  At  pressures  higher  than  0.6  Torr,  the  beam  defocused 
because  of  an  increase  in  beam  current  neutralization  [14]. 

The  above  research  included  Usov’s  well-known  experiment  with  a  1 
MeV  beam  in  a  large  gas -filled  drift  chamber.  In  a  homogeneous 
transverse  magnetic  field,  the  radius  of  curvature  of  the  beam 
trajectory  was  determined  by  the  0.5  MeV  nose  of  the  pulse,  and  not  by 
its  1  MeV  peak  energy,  thus  indicating  that  the  beam  could  track  the 
plasma  channel. 

The  researchers  at  IYaF-TPI  stress  the  importance  of  the  effect  of 
collective  fields  generated  by  the  beam  itself  and  the  beam’s  interaction 
with  the  plasma  channel.  These  researchers  also  found  that  an  IREB 
can  be  effectively  transported  across  a  magnetic  field  along  localized 
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r,  ns 


Fig.  27 — Pulse  length  as  a  function  of  background  pressure  and 
transport  length,  L  [17] 


plasma  channels  having  diameters  smaller  than  that  of  the  electron 
beam  as  well  as  along  thin  conducting  wires.  The  effectiveness  of 
transport  was  found  to  be  determined  by  the  degree  of  magnetic  and 
electrostatic  neutralization  of  the  beam  [14]. 


III.  IAE  RESEARCH  ON  HOSE  INSTABILITY 


Researchers  at  the  Kurchatov  Atomic  Energy  Institute  (IAE),  under 
the  leadership  of  L.  I.  Rudakov,  carried  out  experiments  from  1974  to 
1978  using  the  MS  accelerator  at  0.4  MeV,  40  kA,  40  ns  to  investigate 
beam  transport  and  the  formation  of  plasma  channels  with  the  help  of 
external  magnetic  fields  for  beam  focusing.  The  experiments  demon¬ 
strated  poor  beam  transport  at  background  air  pressures  higher  than 
0.5  Torr  [18,19]. 

Since  1982,  IAE  researchers  have  studied  the  resistive  hose  instabil¬ 
ity  of  an  IREB  beam  injected  into  air  and  gases  [20,21].  Resistive  hose 
instability  is  the  result  of  resonant  interaction  of  magnetic  field  oscilla¬ 
tions  with  betatron  oscillations  of  beam  particles  and  is  observed  as 
macroscopic  transverse  beam  oscillations,  which  eventually  deflect  the 
beam  completely  off  axis  and  onto  the  wall  of  the  drift  tube. 

The  formation  time  of  hose  instability  was  found  to  be  considerably 
shorter  than  the  characteristic  time  of  plasma  magnetic  field  decay. 
The  distance  at  which  the  instability  begins  to  oscillate  is  of  the  order 
of  several  betatron  oscillation  wave  lengths.  The  time  and  length  of 
the  instability  formation  decreased  with  increasing  pressure  from  1  to 
100  Torr. 

Experiments  with  resistive  hose  instability  were  carried  out  using 
the  “Neptun”  accelerator  at  450  keV,  20  kA,  and  pulse  length  of  60  ns. 
The  beam  electrons  passed  through  a  20  micrometer  thick  titanium  foil 
and  were  injected  into  a  drift  chamber  filled  with  N2,  H2,  or  air  at  vari¬ 
ous  pressures.  The  beam  transport  distance  was  varied  from  3  cm  to  3 
m  and  the  pressure  within  the  chamber  was  varied  from  1  to  100  Torr 
[20,21]. 

The  measurements  involved  radial  distribution  of  total  beam  current 
density  at  different  distances  from  the  point  of  beam  injection  at  dif¬ 
ferent  pressures.  Typical  scans  are  shown  in  Fig.  28  [20].  At  1  Torr, 
the  current  density  distribution  maintained  its  axial  symmetry  to  a  dis¬ 
tance  of  0.8  m  (Fig.  28a).  As  beam  transport  distance  increased,  this 
symmetry  was  perturbed  (Fig.  28b)  first  at  the  trailing  edge  of  the 
pulse.  In  some  cases,  beam  deflection  was  less  than  the  beam  radius  at 
a  distance  of  2  meters,  but  at  3  meters  at  the  pulse  trailing  edge,  the 
beam  was  deflected  by  more  than  the  beam  radius.  As  gas  pressure 
increased,  the  asymmetiy  of  the  beam  current  density  profile  and  beam 
deflection  also  increased  (Fig.  28c). 
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a  -  p  =  1  Torr,  L  =  26  cm 
b  -  p  =  1  Torr,  L  =  85  cm 
c  -  p  =  5  Torr,  L  =  85  cm 

Fig.  28 — Beam  current  density  radial  distribution  at  different  times 

after  injection  [20] 


The  measurements  [17,18]  showed  that  the  distance  at  which  the 
beam  is  thrown  from  the  axis  into  the  chamber  wall  depends  upon 
pressure  (Fig.  29)  and  decreases  with  increasing  pressures. 

The  data  obtained  from  Faraday  cups  were  in  good  agreement  with 
optical  observations  of  plasma  scintillation  in  the  drift  chamber 
obtained  from  electron-optical  converter  photographs  (Fig.  30). 

Figure  30a  shows  three  consecutive  photographs  of  the  beam  thrown 
off  axis  onto  the  chamber  wall  at  30  Torr.  In  10  to  15  ns,  the  beam 
remained  close  to  the  axis  of  the  chamber;  in  15  to  20  ns,  the  beam 
moved  radially  through  about  one  radius;  and  in  20  to  25  ns,  the  beam 
was  thrown  completely  onto  the  chamber  wall.  At  60  Torr  (Fig.  30b) 
the  beam  is  thrown  to  the  chamber  wall  in  only  11  ns.  The  charac¬ 
teristic  excitation  length  was  45  cm  at  30  Torr  and  30  cm  at  60  Torr. 
The  average  transverse  displacement  velocity  was  found  to  increase 
with  pressure,  changing  from  5  x  108  to  2  x  109  cm/s,  for  a  pressure 
change  5  to  60  Torr,  respectively  [20]. 

The  experiments  demonstrated  that  the  beam  trajectory  instability 
develops  in  a  spiral  that  unwinds  toward  the  end  of  the  drift  tube  and 
that  beam  asymmetry  is  accompanied  by  beam  rotation  around  its  axis. 
The  angular  velocity  of  this  rotation  was  3  x  10s  sec-1  and  was 
observed  in  both  left  and  right  directions.  The  instability  was  found  to 
possess  a  weakly  convective  character  with  an  oscillation  drift  velocity 
less  than  102  cm/s.  At  pressures  greater  than  10  Torr,  the  transverse 
oscillation  amplitudes  were  found  to  be  so  high  that  a  large  part  of  the 
beam  was  thrown  to  the  chamber  wall. 
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Fig.  29 — Distance  to  off-axis  displacement  as  a  function 
of  pressure  [20] 


a  -  p  =  30  Torr,  L  =  1 .3  m  The  arrow  represents  the  direction  of  the 

b  -  p  -  60  Torr,  L  ^  1.3  m  beam  injection.  The  horizontal  dotted  lines 

represent  the  position  of  the  chamber  wall. 
The  vertical  dark  lines  are  axial  distance 
markers  located  20  cm  apart 

Fig.  30— Electron-optical  converter  graphs  of  axial  displacement  [20] 
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IAE  researchers  state  that  the  accuracy  of  their  time  measurement 
in  the  observations  of  plasma  scintillation  due  to  the  propagation  of 
the  beam  was  within  1  ns.  The  plasma  channel  formed  by  the  beam 
was  observed  to  remain  for  more  than  100  ns  after  the  passage  of  the 
electron  beam  pulse  [20]. 

The  current  density  distribution  in  the  cross  sectional  plane  of  the 
transport  chamber  was  measured  as  a  function  of  time  (Fig.  31)  [20], 
The  figure  shows  the  asymmetry,  displacement,  and  rotation  of  the 
beam  in  the  cross  sectional  view.  Oscillograms  of  signals  from  a  beam 
collector  placed  at  the  location  where  the  beam  was  thrown  from  the 
axis  (at  L  =  85  cm)  at  p  =  15  Torr  are  shown  in  Fig.  32  [20]. 


i=20ns  jgns 


Fig.  31 — Beam  current  density  distribution  (in  A/cm2)  as  a  function  of 

time  [20] 


The  beam  propagated  symmetrically  during  the  first  15  ns;  for  the 
next  10  ns  the  beam  disappeared  from  the  collector’s  field  of  view  and 
reappeared  on  one  side  of  the  collector’s  face.  The  beam  was  displaced 
radially  with  a  velocity  of  109  cm/s. 

The  time  of  plasma  magnetic  field  decay  in  these  pressure  ranges 
was  longer  than  100  ns  for  a  3  cm  beam  radius  and  plasma  conduc- 
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p  =  15  Tort,  L  =  85  cm 

Fig.  32 — Oscillograms  of  the  voltage  V,  diode  current  I,  and  beam 
current  density  Jb  as  measured  with  sectioned  Faraday  cups 
located  at  different  distances  from  the  chamber  axis  [20] 

tivity  tr  >  lO^s-1.1  The  beam  instability  developed  in  considerably  less 
time  (t  =  1  to  20  ns)  than  the  magnetic  field  decay  time. 

Recent  theoretical  work  at  IAE  on  the  resistive  hose  instability  led 
to  the  conclusion  that  under  conditions  of  above  experiments  [20,21], 

1In  the  investigation  of  the  plasma  magnetic  field  decay  time,  Soviet  researchers  do 
not  differentiate  between  the  monopole  time  constant  jm  and  the  dipole  time  constant 
Td.  The  time  constant  to  which  the  Soviets  refer  with  a  value  of  100  ns  is  probably  rm, 
which  is  observed  for  the  plasma  current  at  the  tail  of  a  beam  current  pulse.  The  value 
of  the  dipole  decay  time  constant  rd  would  be  5  or  10  times  shorter.  The  hose  instability 
builds  up  on  the  time  scale  of  rd. 
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an  external  axial  magnetic  field  of  6  kG  applied  to  the  IREB  can 
totally  dampen  short-wave  oscillations  causing  hose  instability.  The 
long-wave  oscillations  do  not  have  time  to  develop  for  the  meter  propa¬ 
gation  length  conditions  in  these  experiments  [22]. 

Theoretical  analysis  was  also  applied  to  beam  self-focusing  produced 
by  two  effects  in  plasma:  (1)  increase  of  beam  space  charge  neutraliza¬ 
tion  with  increasing  plasma  density  and  (2)  generation  of  magnetic 
field  by  thermo-convective  currents  in  inhomogeneous  plasma.  IAE 
researchers  made  a  calculation  to  demonstrate  the  degree  of  compres¬ 
sion  of  a  beam  with  parameters  of  5  MeV,  60  kA,  r  -  3  cm,  in  a  plasma 
with  a  density  of  1016  cm"3  that  extends  along  the  beam  direction  for  a 
distance  of  L  =*  50  cm.  The  generation  of  the  magnetic  field  during 
plasma  heating  in  the  beam  relaxation  process  was  found  to  lead  to 
beam  compression  with  a  compression  factor  of  10  [23]. 

Other  recent  theoretical  work  at  IAE  included  the  research  on  the 
resonant  helical  instability  of  an  IREB  [24]  and  the  study  of  the  gen¬ 
eration  of  a  magnetic  field  in  combination  with  the  filamentation  of  an 
IREB  propagating  through  a  plasma  with  strong  current  neutralization 

[25] .  These  studies  showed  that  for  typical  experimental  parameters 
for  IREB  passing  through  gases  and  plasma  (for  beam  currents  from  30 
to  100  kA,  rb  =  1  to  2  cm,  and  the  plasma  density  of  about  1014)  helical 
instability  is  present  and  fast  oscillations  are  excited  causing  local 
compression  of  the  beam  [25]. 

The  expulsion  of  magnetic  field  from  the  electron  beam  by  hot  elec¬ 
tron  gas  whose  pressure  is  greater  than  the  pressure  of  the  magnetic 
field  has  also  been  recently  investigated  [26].  In  this  theoretical 
analysis,  the  expulsion  process  in  a  high -conductivity  medium  was 
shown  to  be  accompanied  by  a  rearrangement  of  the  magnetic  field  in  a 
time  considerably  shorter  than  skin  effect  time  rsk  {rsk  =  47ror2/c2, 
where  a  =  conductivity,  r  -  characteristic  size  of  the  current  flow  area) 

[26] ,  Nonlinear  wave  solutions  have  been  obtained  for  a  constant 
dispersion  velocity,  which  can  be  considerably  lower  than  the  velocity 
of  fast  electrons. 


IV.  FIAN-IOF  RESEARCH— EXPERIMENTAL 
BEAM  PROPAGATION  STUDIES 


During  the  early  1970s,  experimental  work  at  the  Lebedev  Physics 
Institute  (FIAN)  involved  general  properties  during  beam  transport 
through  gases.  Since  1981  FIAN,  together  with  the  newly  formed 
Institute  of  General  Physics  (IOF)  under  the  leadership  of  A.  A. 
Rukhadze,  has  been  publishing  experimental  studies  of  beam  current 
enhancement,  thermal  effects  within  the  IREB,  and  their  effect  on  the 
background  gas  density. 

The  first  experimental  results  concerning  the  perturbation  of  air 
density  by  an  IREB  pulse  were  published  by  FIAN-IOF  in  1984  [27]. 
The  gas  dispersion  after  the  initiation  of  the  beam  pulse  was  analyzed 
numerically  [28].  In  the  latest  experiment  beam  parameters  were 
measured  and  gas  density  change  was  determined  along  the  beam  axis 
[29].  A  simple  relationship  was  found  between  the  measured  gas  den¬ 
sity  distribution  after  the  beam  pulse  and  the  input  energy  profile  [29], 

The  behavior  of  gas  density  after  the  injection  of  an  IREB  pulse  was 
observed  [27]  experimentally  using  the  Terek-IR  accelerator  [30,31] 
with  beam  currents  of  10  to  15  kA,  beam  energies  of  1.5  MeV,  and 
pulse  lengths  of  about  60  ns.  The  current  densities  on  axis  were 
0.3  kA/cm2  and  1  kA/cm2  with  the  beam  radius  of  1.5  to  3  cm.  The 
beam  was  injected  into  atmospheric  air  through  a  50  pm  thick  titanium 
foil. 

The  gas  density  radial  distribution  for  different  times  after  beam 
injection  is  represented  in  Fig.  33  .[27].  There  were  considerable  differ¬ 
ences  in  the  behavior  of  background  gas  with  different  beam  current 
density  values.  This  type  of  gas  behavior  could  be  explained  simply  by 
the  difference  in  energy  input;  FIAN-IOF  researchers  attribute  it  to  the 
energy  accumulated  in  the  internal  degrees  of  freedom  of  the  gas  and 
to  the  strong  dependence  of  relaxation  time  on  gas  temperature.  The 
relaxation  time  was  said  to  be  determined  by  the  input  energy  of  the 
beam.  This  gas  behavior  can  also  be  seen  in  the  dependence  on  time  of 
the  rarified  gas  density  along  the  beam  axis  after  beam  injection  as 
shown  in  Fig.  34  127]. 

The  dependence  of  the  maximum  amplitude  of  the  gas  compression 
coordinate  upon  the  time  after  the  beam  injection  with  a  current  den¬ 
sity  along  the  beam  axis  of  j  =1  kA/cm2  is  represented  in  Fig.  35  127]. 

These  experiments  determined  that  the  velocity  of  the  compression 
wave  outside  the  volume  occupied  by  the  beam  itself  (r  >:  3  cm)  is  close 
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2  -  j  =  1  kA/cm2 

Fig.  34 — Gas  density  along  the  axis  of  the  beam  as  a  function  of  time 
after  the  IREB  injection  [27] 


tjiiS 

Fig.  35— Maximum  amplitude  of  the  gas  compression  coordinate  as  a 
function  of  time  after  the  beam  injection  for  j  =  1  kA/cm2  [27] 
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to  the  velocity  of  sound.  The  investigation  of  the  density  behavior 
after  IREB  injection  can  be  a  valuable  source  of  information  about  the 
relaxation  kinetics  of  the  oscillation-induced  conditions. 

The  decrease  in  gas  density  along  the  beam  trajectory  was  further 
experimentally  investigated  [29]  using  the  Terek-IR  accelerator  with 
beam  currents  from  4  to  15  kA  at  beam  energies  from  1  to  1.2  MeV 
and  pulse  lengths  of  about  60  ns.  A  relationship  was  determined 
between  gas  density  distribution  after  the  IREB  pulse  and  the  energy 
deposition  profile.  A  typical  gas  density  distribution  5p(r)  for 
t observed  =  34  /is  is  shown  in  Fig.  36a  [29].  The  radial  distribution  of 
the  beam  energy  release  into  the  gas  is  shown  in  Fig.  36b  [29]. 


Solid  line  -  calculation 

Data  Points  •  +  -  experimental  observation 


Fig.  36 — a  -  Gas  density  radial  distribution  observed  34  /is  after  the 
IREB  pulse  [29] 

b  -  Radial  distribution  of  the  beam  energy  release  into  the  gas 


The  gas  density  change  along  the  beam  axis  as  a  function  of  time  is 
shown  in  Fig.  37  [29].  The  minimum  value  of  gas  density  near  the 
beam  axis  can  be  seen  from  Fig.  38  where  its  dependence  on  p/p0  is  in 
good  agreement  with  theoretical  calculations.  Another  issue  was  the 
total  beam  current  enhancement  when  the  plasma  electron  current 
flows  in  the  same  direction  as  the  IREB  and  its  total  amplitude  can 
exceed  the  amplitude  of  the  injected  beam.  Researchers  at  Yefremov 
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Solid  line  -  theory 
o  -  experimental  observations 

Fig.  37 — Dependence  of  the  gas  density  along  the  beam  axis 
upon  time  [29] 


Solid  line  represents  calculated  values 

Data  points  represent  experimental  observations  made  in  Ref.  27 

Fig.  38 — Dependence  of  the  released  energy  density  along  the  beam 
axis  upon  the  change  in  gas  density  at  tobgerved  »  t8  [29] 
here:  ts  ^  where  a  =  effective  beam  radius 
C$  -  velocity  of  sound 
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(NIIEFA)  observed  the  amplification  of  the  beam  current  in  1978-1979 
[32,33]  and  they  were  able  to  correlate  this  experimental  finding  with 
the  movement  of  plasma  electrons  in  the  turbulent  fields  of  a  plasma- 
beam  instability. 

The  current  enhancement  was  experimentally  determined 
[32,33,34,35]  to  be  due  to  two-stream  instability.  However,  in  1984, 
FIAN-IOF  researchers  observed  beam  current  enhancement  at  pres¬ 
sures  close  to  atmospheric,  contrary  to  what  these  experiments  led 
them  to  believe  [36]. 

Again  using  the  Terek-lR,  with  7  kA,  1  MeV,  and  60  ns  pulse 
length,  the  beam  was  injected  into  air,  He,  SF6,  and  a  mixture  of 
SF6  +  He  +  H2.  The  experiments  were  carried  out  both  with  an  axial 
magnetic  field  (B  -  2.4  kG)  and  without  a  magnetic  field. 

The  experiments  have  shown  nanosecond  spikes  with  amplitudes 
several  times  higher  than  that  of  the  injected  beam.  This  effect  was 
observed  for  all  gases  at  all  distances  from  the  point  of  injection. 

The  characteristic  length  of  the  spike  in  the  total  beam  current 
pulse  was  found  to  correspond  to  the  time  of  beam  neutralization  [30]. 
This  was  assumed  to  be  the  effect  of  current  enhancement  at  the  beam 
front.  The  spike  structure  of  the  pulse  was  visible  with  pressure  start¬ 
ing  from  several  tens  of  Torr  at  distances  greater  than  1  meter  with 
the  spike  amplitude  increasing  toward  the  end  of  the  pulse  [36],  The 
application  of  the  axial  magnetic  field  completely  obliterated  this  effect 
even  at  pressures  close  to  atmospheric.  The  effect  was  most  prominent 
in  helium,  but  it  was  also  observed  in  the  other  gases. 

The  effect  of  total  beam  current  enhancement  was  attributed  to  the 
appearance  of  the  spikes  in  the  current  pulse.  At  fairly  high  pressures 
and  at  transport  distances  of  greater  than  1  meter,  at  the  end  of  the 
total  current  pulse  a  nanosecond  peak  was  observed  that  was  a  factor 
of  1.5  greater  than  the  amplitude  of  the  injected  beam  current.  This 
effect  of  beam  current  enhancement  increased  with  transport  distance 
and  with  increase  in  background  gas  pressure  (Fig.  39)  [36], 

In  He,  even  at  p  =  750  Torr,  the  peaks  could  not  be  dampened  with 
the  axial  magnetic  field,  although  other  gases  did  not  show  this  effect. 
There  was  thus  a  difference  in  the  conductance  of  the  beam  plasma  in 
He  than  in  the  other  gases.  FIAN-IOF  researchers  are  at  present 
investigating  this  beam  current  enhancement  theoretically  [36]. 

Other  experiments  performed  with  Terek- 1R  included  plasma  return 
current  induced  by  IREB  and  concentration  and  temperature  of  beam- 
produced  plasma,  using  laser  dispersion  [37].  At  pressures  greater  than 
3  Torr,  plasma  concentration  increased  to  5  x  1016  cm-3  at  a  tempera¬ 
ture  of  T  ~  1  eV,  which  demonstrates  the  presence  of  collective 
processes  that  take  place  in  the  plasma.  For  these  experiments  the 
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1  l  I.  ,1  .1  I 
50  ns 


Solid  line  B0  «  0 
Dotted  line  B0  =  2.4  kG 


Fig.  39— a  -  Total  beam  current  as  a  function  of  distance  from 
injection  at  a  pressure  of  75  Torr  in  He  [36] 
b  -  Total  beam  current  as  a  function  of  pressure  in  He  at  a 
distance  of  113  cm  from  injection 


accelerator  was  operated  with  a  beam  current  of  25  kA  at  1.5  MeV  with 
a  60  ns  pulse  length. 

The  Terek-lR  was  used  in  1984  in  a  set  of  experiments  on  the 
investigation  of  beam  pulse  front  erosion  by  neutral  gases  at  pressures 
from  100  to  750  Torr  [38],  These  experiments  were  performed  at  1.5 
MeV,  12  kA,  40  ns  using  a  drift  tube  1.2  m  long.  The  gases  used  were 
nitrogen  and  a  mixture  of  N2  and  SF6  (in  the  ratio  of  6  to  1.)  The 
beam  transport  through  the  drift  tube  was  carried  out  with  and  without 
the  application  of  an  axial  magnetic  field  of  about  5.8  kG.  Erosion  of 
the  leading  front  of  the  beam  pulse  was  observed  to  increase  with  beam 
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transport  length.  At  a  distance  of  60  cm,  the  pulse  width  was  short¬ 
ened  to  half  its  original  width.  With  the  application  of  the  axial  mag¬ 
netic  field  this  effect  of  pulse  front  erosion  completely  disappeared. 
The  erosion  of  the  pulse  front  was  observed  over  the  whole  pressure 
range  from  100  to  750  Torr,  and  any  quantitative  effect  of  the  erosion 
was  found  not  to  depend  upon  gas  pressure  [38], 

Another  set  of  experiments  was  performed  in  1984  on  the  injection 
of  IREBs  into  gas  at  these  high  pressures  but  with  lower  beam  current 
densities  (/&  <  50  A  A cm2).  The  Terek-2  accelerator  was  used  for  these 
investigations  with  beam  energy  of  360  keV  and  beam  currents  of  1.4 
kA  and  250  A  [39], 

In  experiments  performed  earlier  at  FIAN  using  the  Terek-lR  [30], 
a  dependence  of  the  plasma  return  current  upon  the  beam  current  as  a 
function  of  background  gas  pressures  was  obtained.  In  subsequent 
experiments  [37]  the  reverse  current  was  found  to  have  a  maximum  in 
the  neighborhood  of  p  ~  1  Torr  just  where  a  considerable  concentra¬ 
tion  of  plasma  electrons  was  found  (Fig.  40)  [37].  As  pressure 
increased,  the  role  of  cumulative  ionization  (ionization  by  secondary 
electrons)  decreased  and  direct  gas  ionization  by  the  IREB  became 
responsible  for  the  formation  of  the  plasma. 

The  FIAN-IOF  researchers  stressed  the  significance  of  their  research 
to  optimal  IREB  transport  through  gas,  as  well  as  to  high-power  gas 
lasers  and  plasmo-chemical  reactors.  The  results  were  also  claimed 
applicable  to  IREB  diagnostics  in  gas,  where  the  measurements  of  the 
electron  beam  parameters  give  no  single-valued  results  or  perturb  the 
distribution  of  induced  currents  or  fields. 


+  -  He 
•  -air 
o  -SF6 

Fig.  40 — Plasma  current  to  beam  current  ratio  as  a  function  of  gas 

pressure  [37] 


V.  KhFTI  EXPERIMENTAL  RESEARCH  ON 
BEAM  PROPAGATION 


KhFTI’s  experimental  research  on  the  propagation  of  IREBs 
through  background  gas  or  plasma  can  be  broken  down  into  two 
separate  groups.  One  group,  whose  publications  extend  from  1973  to 
1982  under  the  leadership  of  Yu.  V.  Tkach  and  Ya.  B.  Faynberg,  has 
investigated  IREB  propagation  through  air  using  a  1  MeV,  50  kA,  30 
ns  electron  beam  with  transport  distances  of  up  to  3  m.  The  second 
group,  headed  by  Ya.  B.  Faynberg  with  A.  K.  Berezin  as  a  leading 
author,  and  publishing  from  the  early  1970s  to  1981,  used  a  2  MeV  and 
a  20  MeV  linear  electron  accelerator  with  low  beam  currents  of  0.5  to 
1.0  A  and  investigated  the  interaction  of  IREB  with  air  and  plasma. 
No  papers  of  KhFTI  researchers  have  been  found  since  1982. 

The  early  work  of  the  first  group,  dating  back  to  1973,  dealt  with  the 
formation  of  return  current,  the  optimum  pressure  for  transport  in  air, 
and  the  measurement  of  the  ionization  front  velocity  in  various  gases 
in  the  pressure  range  from  5  x  10-2  to  10  Torr  [40].  In  1980,  the 
effect  of  collisions  upon  the  beam  relaxation  and  of  collective  and 
elementary  processes  upon  energy  transport  through  gas  were  studied 
[41].  The  later  experiments  used  a  0.7  MeV,  30  to  35  kA  electron 
beam,  40  ns  pulse  length,  with  transport  distances  up  to  3  m  in  air,  in 
the  pressure  range  from  0.1  to  100  Torr,  without  the  presence  of  exter¬ 
nal  magnetic  fields. 

The  KhFTI  researchers  claim  to  have  demonstrated  experimentally 
that  the  relaxation  of  an  IREB  in  a  neutral  gas  and  plasma  without 
external  magnetic  fields  occurs  mainly  through  collective  processes, 
which,  during  their  formation,  are  accompanied  by  the  excitation  of  an 
intensive  microwave  and  X-ray  radiation  from  the  plasma,  the  modula¬ 
tion  of  the  IREB,  and  the  formation  of  a  fast  group  of  particles  with 
energies  reaching  20  to  30  keV.  According  to  computer  simulation,  in 
addition  to  the  collective  processes,  an  important  role  is  played  by  elas¬ 
tic  scattering  of  the  beam  particles  on  the  plasma  atoms  and  ions, 
which  removes  particles  from  the  beam  channel  if  the  scattering  angle 
is  greater  than  a  certain  critical  value. 

The  KhFTI  researchers  devised  complex  diagnostics  for  these  exper¬ 
iments  to  obtain  accurate  measurements.  A  series  of  Rogowski  loops 
monitored  the  diode  and  beam  currents  at  various  locations  along  the 
drift  tube,  a  25-strip  calorimeter  measured  the  total  beam  energy  and 
its  radial  distribution  and  allowed  control  of  the  beam  diameter  while 
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making  calorimetric  measurements,  a  12-channel  magnetic  analyzer 
determined  the  spectral  distribution  of  the  beam  electrons,  and  a  4- 
channel  X-ray  spectrometer  measured  bremsstrahlung  radiation 
emanating  from  the  plasma  heated  by  the  IREB. 

The  KhFTI  researchers  stressed  the  importance  of  transporting 
energy  to  large  distances  and  mentioned  that  the  most  important 
characteristic  is  the  beam  energy  density  obtained  on  target  as  a  func¬ 
tion  of  background  gas.  The  measured  energy  carried  by  the  beam  over 
different  distances  as  a  function  of  pressure  is  shown  in  Fig.  41..  The 
energy  carried  by  the  beam  as  a  function  of  distance  for  various  pres¬ 
sures  is  shown  in  Fig.  42..  The  condition  for  maximum  energy  carried 
by  the  beam  to  a  certain  distance  is  shown  in  Fig.  43  as  a  function  of 
pressure  [41]. 

From  their  experimental  measurements,  the  Soviet  researchers  con¬ 
cluded  that  in  order  to  obtain  the  maximum  energy  carried  by  the 
beam  for  short  distances  (  <  50  cm)  it  is  necessary  to  choose  the  pres¬ 
sure  at  which  the  beam  is  self-focusing  at  p  «  1  Torr.  To  carry  the 
beam  over  long  distances,  it  is  necessary  to  take  beam  energy  dissipa¬ 
tion  into  account.  The  optimum  gas  pressures  for  long  distance  beam 
transport  were  found  between  0.2  and  0.6  Torr  [41]. 

The  changes  in  particle  velocity  and  particle  density  in  beam  propa¬ 
gation  through  gas  were  measured  experimentally.  The  beam  energy 
spread  was  measured  at  the  output  of  the  accelerator  and  at  various 
locations  along  the  drift  tube  as  a  function  of  gas  pressure  and  was 
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lb  =  30  kA,  rb  =  1 .8  cm,  Rc  =  5.5  cm;  z  is  path  length: 

1  -  z  =  6  cm  3  -  z  =  1  m 

2  -  z  =  .5  m  4  -  z  =  3  m 

eq  =  maximum  measured  beam  energy 

Fig.  41 — Energy  carried  by  the  beam  to  different  distances  as  a 
function  of  gas  pressure  [41] 
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lb  =  30  kA,  rb  =  1 .8  cm,  =  5.5  cm 

1  -  p  =  0.1  Torr  3  -  p  =  1  Torr 

2  -  p  =  0.4  Torr  4  -  p  =  10  Torr 


Fig.  42 — Energy  carried  by  the  beam  as  a  function  of  distance  for 
various  gas  pressures  [41] 


lb  =  30  kA,  rco|iector  =  3  cm,  rb  =  1.8  cm, 
Rc  =  5.5  cm 


Fig.  43 — Pressure  required  for  maximum  energy  carried  by  the  beam  to 
different  distances  [41] 
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represented  as  a  time  function  for  two  pressures  in  Fig.  44  [41].  At  the 
optimum  transport  pressure  of  p  =  0.4  Torr,  beam  energy  distribution 
changed  only  slightly  with  change  in  transport  distance  (see  Fig.  44). 
At  high  pressures,  the  beam  energy  peak  was  shifted  to  the  lower  ener¬ 
gies  and  the  low  energy  particles  decreased  with  increased  z,  while  the 
pulse  shortened. 

The  beam  energy  spread  was  determined  as  a  function  of  gas  pres¬ 
sure  for  several  transport  distances  (Fig.  45)  [41].  The  energy  spread 
was  independent  of  p  and  z  at  short  transport  distances,  but  it  was 
markedly  dependent  at  the  3  m  point. 

Further  analysis  of  their  data  allowed  the  researchers  to  calculate 
the  specific  beam  particle  energy  as  shown  in  Fig.  46  demonstrating  a 
much  larger  beam  energy  loss  in  the  pressure  region  of  0.1  and  1  Torr 
than  in  the  region  of  0.3  to  0.6  Torr. 

The  researchers  noted  that  at  the  pressure  range  of  0.1  to  0.2  Torr, 
the  increase  in  energy  loss  was  due  to  incomplete  space  charge  neutral¬ 
ization  and  deceleration  of  the  beam  by  the  space  charge  field.  The 
increased  loss  in  the  1  to  3  Torr  range  was  due  to  collective  processes 
accompanied  by  microwave  radiation.  The  dependence  of  microwave 


p  =  0.4  Torr 


t,  ns 


n,  rel.  units  tr  ns 


p  =  4  Torr 


z  =  6  cm 


z  =  0.3  m 


z  =  3  m 


BO  210310  UZO  550  B^keV 


Fig.  44 — Time  plots  showing  the  IREB  energy  spread  at  three  locations 
along  the  drift  tube  for  p  =  0.4  and  4  Torr  [41] 
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n,  reL  units 


1  -  p  =  0.1  Torr  2  -  p  =  0.4  Torr 

3  -  p  =  1  Torr 


Fig.  45 — Beam  energy  spectra  as  a  function  of  background  pressure  at 
two  transport  distancesjjtlj 


Ey/n,  reL  units 


1  -  z  =  6  cm  2  -  z  =  30  cm  3  -  z  =  3  m 

here  Eb  is  beam  energy  and  n  is  the  number  of  beam  particles 


Fig.  46 — Specific  energy  of  the  transported  electron  beam  at  different 
transport  distances  as  a  function  of  gas  pressure  [41] 
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radiation  intensity  upon  pressure  and  transport  length  is  shown  in  Fig. 
47  for  the  wavelengths  of  3  cm,  8  mm,  and  1.8  mm. 

Figure  47  shows  that  a  microwave  radiation  peak  appears  at  points 
in  the  pressure  range  where  the  beam  sustains  a  considerable  particle 
energy  loss.  The  microwave  radiation  shows  a  minimum  at  the  point 
of  optimum  beam  transport.  As  the  beam  path  lengthens,  the  range  of 
pressures  for  which  microwave  radiation  is  observed  increases. 


microwave  radiation  wavelength: 

-  .  -  4  -  A  -  3  cm 

.  . .  \  =  3  mm 

—  —  —  — - A  -1.8  mm 

Fig.  47 — Dependence  of  the  microwave  radiation  intensity  upon  gas 
pressure  and  transport  length  [41] 


X-ray  radiation  emanating  from  the  plasma  was  found  to  be  fairly 
weak  in  the  pressure  range  corresponding  to  optimum  beam  transport 
and  to  increase  at  pressures  from  0.6  to  3  Torr  in  the  region  where  the 
maximum  modulation  of  the  beam  current  was  observed.  The  depen¬ 
dence  of  the  X-ray  radiation  intensity  as  a  function  of  gas  pressure  is 
shown  in  Fig.  48  [41]. 

The  angular  and  radial  beam  current  distribution  as  a  function  of 
gas  pressure  and  time  were  investigated  [42]  using  a  beam  current  of 
30  kA,  a  beam  energy  of  0.5  MeV,  pulse  length  of  45  ns,  beam  radius  of 
1.8  cm,  and  background  gas  pressure  from  5  x  10“2  to  50  Torr.  The 
findings  are  shown  in  Fig.  49.  The  initially  cylindrical  radial  electron 
distribution  assumed  a  hollow-beam  character  at  a  distance  of  40  cm 
from  the  point  of  beam  injection.  At  low  pressures  (less  than  2  Torr) 
this  happens  in  a  time  equivalent  to  the  pulse  length;  at  higher  pres- 
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lR,  rel.  units 


Fig.  48— Dependence  of  the  soft  X-ray  radiation  (15  to  35  keV) 
emanating  from  the  plasma  upon  background  gas  [41] 


a,  degrees 


1  -  15  ns,  2  -  30  ns,  3  -  50  ns 
I  -  7  x  1 0‘2  Tor r  II  -  2  x  10’1  Torr  III  -  2  Torr 


Fig.  49 — Radial  distribution  of  beam  current  density  (a)  and  IREB 
particle  angle  distribution  (b)  in  the  drift  tube  for  various 
times  and  pressures  [42] 
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sures,  the  beam  possesses  a  well-defined  hollow  beam  structure  during 
the  entire  pulse. 

Using  the  2  MeV  and  the  20  MeV  linear  accelerators  with  a  current 
range  of  0.5  to  1.0  A,  the  second  group  of  KhFTI  researchers  studied 
the  interaction  of  an  IREB  with  plasma.  The  early  experimental 
research  of  this  group  (1971-1972)  concentrated  on  IREB  interaction 
with  low  density  hydrogen  plasma  of  1010  to  1013  cm-3  using  the  2  MeV 
electron  linear  accelerator  at  low  beam  currents  of  about  1  A 
[43,44,45].  Later  experiments  (1974  to  1976)  involved  high  density 
plasma  of  1016  to  1017  cm-3  [46,47,48].  In  experiments  carried  out  in 
1979  with  the  2  MeV  IREB  in  air  at  atmospheric  pressure,  KhFTI 
researchers  demonstrated  the  formation  of  a  intermediately  dense 
plasma  (up  to  1013  cm”3)  and  also  observed  an  effective  collective 
interaction  of  the  electron  beam  with  plasma  [49].  The  latter  was 
apparent  from  the  considerable  loss  of  beam  energy  as  well  as  from 
microwave  and  soft  X-ray  radiation  localized  in  the  area  of  beam 
interaction. 

In  1978-1979  the  group  published  two  papers  on  the  LU-40  linear 
accelerator  with  20  MeV,  0.5  A,  and  pulse  length  of  10  ns  [50,51]. 
Only  about  1  percent  of  the  beam  passed  through  plasma  with  a  den¬ 
sity  of  6  to  8  x  1016  cm-3.  At  2  to  3  x  1016  cm-3,  about  20  percent  of 
the  beam  could  be  propagated.  However,  at  a  density  of  5  x  1015  cm-3 
only  about  15  to  20  percent  of  the  beam  would  be  lost. 

The  energy  of  the  X-ray  radiation  emanating  from  plasma  varied 
from  40  to  90  keV,  depending  on  plasma  density,  with  an  intensity  of 
105i?/cm2-s  at  10  cm  from  the  system  axis.  The  Soviet  researchers 
claim  that  the  X-ray  radiation  measurements  correlated  well  with  the 
observed  beam  current  losses. 

The  interaction  zone  between  the  plasma  and  beam  was  less  than  10 
cm;  the  beam  electrons  accelerated  by  an  additional  4  MeV  (for  an  initial 
beam  of  20  MeV).  The  researchers  therefore  proposed  that  high- 
frequency  fields  of  at  least  400  kV/cm  have  been  set  up  in  the  interaction 
region.  They  have  also  deduced  that  the  effectiveness  of  the  collective 
interactions  of  a  monoenergetic,  low  angular  divergence  IREB  increases 
with  the  plasma  density  [50,51]. 


VI.  NIIEFA  EXPERIMENTAL  RESEARCH— 
TWO-STREAM  INSTABILITY 


NIIEFA’s  experiments,  under  the  leadership  of  L.  V.  Dubovoy,  pro¬ 
vided  important  information  on  “reverse  current  enhancement”  (two- 
stream  instability)  in  the  propagation  of  an  IREB  through  air.  The 
research  was  published  in  1977  [32,33]  with  references  to  earlier 
research  dating  back  to  1973  [52,53]  with  no  new  publications  since 
1977. 

A  1.5  MeV,  50  kA,  16  ns  electron  beam  from  the  REP-5  accelerator 
was  injected  into  air  and  nitrogen  at  pressures  from  10-2  to  760  Torr. 
Most  of  the  investigations  were  carried  out  well  below  the  pressure  of  1 
Torr  where  the  two-stream  instability  comes  into  play.  But  measure¬ 
ments  were  also  performed  at  atmospheric  pressure  in  air.  The  charac¬ 
teristics  of  the  propagating  beam  were  measured  at  z  -  60  cm  and 
1.2  m. 

The  investigation  focused  on  the  case  in  which  the  amplitude  of  the 
current  induced  in  plasma  by  the  IREB  pulse  was  higher  than  the 
beam  current  pulse  amplitude.  This  condition  was  called  reverse 
current  enhancement  (or  RCE  effect).  The  RCE  effect  was  present 
only  in  the  case  where  the  electron  beam  has  a  large  energy  spread. 

A  typical  electron  beam  current  pulse  scan  is  shown  in  Fig.  50 
demonstrating  the  beam  current  /*,(£),  total  current  7r(t),  and  reverse 
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Fig.  50 — Beam  current  pulse  scan  [33] 
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current  Ir(t),  where  Ir(t)  =  IT(t)  -  Ib(t)  for  a  pressure  of  0.05  Torr 
[33].  Figure  51  demonstrates  the  dependence  of  the  basic  parameters 
observed  in  this  reverse  current  process  upon  the  background  gas  pres¬ 
sure  [33]. 

From  the  experimental  results,  NIIEFA  researchers  noted  that  the 
reverse  current  Ir  can  be  larger  than  the  beam  current  Ib  with  the 
maximum  (as  shown  for  maximum  y  in  Fig.  51)  found  at  a  pressure  of 
0.03  Torr  marked  by  a  maximum  in  the  beam  energy  spread  (A E)  and 
a  maximum  loss  in  the  energy  of  the  beam. 


(a) 


(b) 


(c) 


•  -  z  =  60  cm  V  -  z  =  1 .2  m 
a  -  beam  energy  E(f  )p,  here  o  -  beam  radius  rb(f)p 
b-Y(p)  =  (lr/lb)(p) 

c  -  AE(p)  change  in  energy  of  the  beam  spectra 

Fig.  51 — Dependence  of  system  parameters  upon  background 
gas  pressure  [33] 


VII.  NIITP  EXPERIMENTAL  RESEARCH— IREB 
INJECTION  INTO  ATMOSPHERIC  AIR 


From  1977  to  1981  the  NIITP  carried  on  a  seemingly  minor  effort 
(represented  by  two  published  papers  and  four  authors)  under  the 
leadership  of  V.  M.  Iyevlev.  A  500  keV  CW  electron  beam  was  injected 
into  air  at  atmospheric  pressure  through  a  differential  pumping  system 
instead  of  a  titanium  foil  window  [54,55],  Special  attention  was  given 
to  the  development  of  the  differential  pumping  system  utilizing  a  series 
of  apertures  with  various  pumping  schemes  incorporated  between  each 
aperture  in  order  to  reduce  the  gas  load  and  provide  a  smooth  transi¬ 
tion  between  the  atmospheric  air  and  the  high  vacuum  required  in  the 
accelerator  section. 

The  maximum  propagation  distance  of  200  keV  electrons  through 
atmospheric  air  was  only  about  30  cm.  However,  as  the  air  was  heated 
up  by  the  beam,  a  lower  density  channel  formed,  extending  the  propa¬ 
gation  distance.  No  further  information  was  published  on  the  details 
of  beam  propagation  through  atmospheric  air. 

Stated  practical  applications  of  the  research  were  beam  welding,  cut¬ 
ting,  and  hole  boring  of  solid  materials. 
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VIII.  VEI  RESEARCH— THEORY  AND 
EXPERIMENT 


VEI  researchers  have  recently  performed  extensive  research  on 
IREB  transport  through  low-pressure  background  gases  (10-3  to  10-5 
Torr),  as  described  in  a  previous  report  [1].  In  high  gas  pressure, 
IREB  propagation  studies  that  began  in  1979  are  still  in  progress.  VEI 
researchers,  under  the  leadership  of  A.  V.  Zharinov  and  M.  A.  Vlasov, 
concentrated  on  the  study  of  the  propagation  windows  for  IREB  trans¬ 
port,  the  determination  of  the  upper  limit  of  gas  pressure  to  maintain  a 
plasma-beam  discharge  (plasma  formed  by  IREB  injected  into  air), 
injection  methods  of  an  IREB  into  the  atmosphere  using  differential 
pumping,  the  formation  of  large  microwave  signals  inherent  in  the 
beam,  and  the  formation  of  positively  charged  electron  beams.  The 
latter  of  these  seems  to  be  of  interest  because  the  researchers  claim 
that  the  beam  does  not  become  subjected  to  any  irregular  instabilities 
and  is  axially  self-centering. 

The  upper  limit  of  the  gas  pressure  for  a  plasma -beam  discharge  was 
experimentally  determined  in  [56]  at  VEI  using  a  300  keV,  20  A,  400 
fis ,  0.6  cm  diameter  electron  beam  and  injecting  the  beam  into  air 
through  a  differential  pumping  system.  The  pressure  range  used  was 
0.1  to  300  Torr.  It  was  determined  that  in  the  pressure  range  50  to  300 
Torr  the  electron  beam  propagation  through  the  dense  air  was  stable 
and  that  the  dynamics  of  the  beam  radius  change  along  the  length  and 
time  of  propagation  can  be  well  approximated  by  the  gas  collision 
theory. 

Researchers  at  VEI  have  investigated  the  decrease  in  the  number  of 
neutral  gas  atoms  in  the  path  of  the  plasma  channel  due  to  the  effect 
of  the  hydrodynamic  expulsion  of  the  channel  hot  gas  [57,56].  Here 
the  interaction  of  the  IREB  with  the  background  gas  has  been  investi¬ 
gated  showing  the  basic  mechanism  to  be  the  dissipation  of  energy  due 
to  electron  collisions  with  the  background  gas  particles  whereby  plasma 
is  formed  and  the  gas  is  heated  up. 

The  gas  temperature  was  estimated  from  the  following  expression 

T  ~  ( Pah )°'2 

where  p0  =  the  initial  gas  density 
lb  =  beam  current. 
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For  pressures  of  10  to  40  Torr,  the  value  of  T  was  found  to  be  0.4  to 
0.5  eV  for  the  observed  beam  parameters;  the  concentration  of  neutrals 
inside  the  channel  was  calculated  to  be  2.5  to  8  x  1016  cm-3  [56]. 

VEI  researchers  carried  out  experiments  in  1984  [58]  on  the  forma¬ 
tion  of  a  beam-plasma  discharge  in  the  area  where  the  electron  beam 
was  injected  into  the  atmosphere  through  a  differential  pumping  sys¬ 
tem.  The  energy  of  the  beam  was  300  keV,  beam  current  27  A,  pulse 
length  500  fis,  and  diameter  0.6  cm.  The  beam  transport  was  stable 
through  the  differentially  pumped  beam  injection  area  where  the 
plasma  density  had  an  extended  gradient  and  the  gas  pressure  changed 
from  20  to  400  Torr. 

The  interaction  of  the  electron  beam  with  the  background  gas  in 
this  area  led  to  the  heating  of  the  gas  and  the  gas  radial  expansion, 
which  generated  a  cylindrical  shock  wave.  The  latter  was  reflected  off 
the  walls,  collapsed  onto  the  axis,  and  radially  expanded.  Because  of 
the  gas  heating,  a  gas  plug  was  formed  in  the  channel  and  obstructed 
the  gas  flow.  It  was  demonstrated  that  because  of  the  conditions  set 
up  by  the  gas  density  gradient  along  the  beam  axis,  a  plasma-beam 
discharge  was  created  with  microwave  radiation  in  the  3  cm  wavelength 
region. 

VEI  researchers  in  1984  also  described  a  differential  pumping  system 
they  have  designed  for  the  injection  of  an  electron  beam  from  a  pres¬ 
sure  of  10-5  Torr  to  a  region  under  pressure  of  10  to  100  Torr  [59]. 

An  interesting  finding  of  the  recent  VEI  research  was  the  formation 
of  overall  positively  charged  beams  of  electrons  [60,61].  These  beams 
had  an  overabundance  of  positive  charge  in  the  neighborhood  of  103  to 
104  V  and  could  be  varied  as  required  by  changing  either  the  pressure 
or  the  type  of  background  gas  or  metal  vapor.  The  critical  current  was 
shown  to  increase  several  hundred  fold  (over  10  A)  using  gridded  con¬ 
ducting  tubes  for  propagation.  The  electric  fields  in  this  case  could  be 
higher  than  in  the  quasi-neutral  case  (when  ni  -  ne9  Ti  =  re)  by  factors 
of  100  to  1000.  The  main  advantage  of  positively  charged  electron 
beams  is  that  they  are  not  subject  to  irregular  instabilities  and  have  no 
beam  current  limits  because  of  these  instabilities  and  are  also  axially 
self-centering,  which  may  allow  them  to  be  easily  transported  along 
various  configurations  of  metal  beam  pipes  [60].  The  formation  of  this 
type  of  electron  beam  is  explained  by  VEI  researchers  using  the 
apparatus  shown  in  Fig.  52  [60].  The  electron  beam  flows  with  a  beam 
current  /*,,  velocity  i >  =  0c  along  the  axis  of  a  tube  of  length  L  and 
radius  R  as  shown  in  the  figure.  A  potential  Vt  is  applied  to  the  tube, 
which  is  negative  with  respect  to  other  electrodes  that  contact  the 
plasma  created  by  the  beam.  All  the  secondary  electrons  created  by 
the  beam  inside  the  tube  have  to  leave  by  the  ends,  and  all  the  ions 
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Fig.  52 — Apparatus  for  the  formation  of  positively  charged  electron 
beams  showing  the  voltage  radial  distribution  [60] 
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must  flow  to  the  walls  under  the  action  of  the  field  Er  where  they 
recombine  at  the  wall  surface.  For  the  secondary  electrons  to  leave  the 
tube,  the  voltage  on  the  tube  VT  must  be  greater  than  the  voltage 
between  the  beam  boundary  and  the  inside  of  the  tube  <pr  (as  shown  in 
the  voltage  distribution  plot  in  Fig.  52).  The  condition  for  the  electron 
beam  to  pass  through  the  tube  must  be  met:  Here  the  kinetic  energy 
of  the  beam  must  be  larger  than  Vy  —  <pr.  The  calculations  in  Ref. 
60  demonstrate  that  this  sort  of  overcompensated  electron  beam  can 
indeed  be  formed  and  could  be  of  future  interest  because  of  its 
aforementioned  advantages. 

VEI  also  considered  the  practical  aspects  of  the  electron  beam  for 
processing  of  materials  [57,62],  involving  the  injection  of  0.5  MeV 
beam  into  air  at  atmospheric  pressure.  It  was  found  that  electron 
beam  propagation  is  very  weakly  dependent  upon  beam  current 
(L  -  1^)  but  strongly  dependent  upon  beam  energy  (L  ~  Efy.  Beam 
energies  of  the  order  of  0.5  MeV  were  recommended  for  processing  of 
materials.  Even  at  this  energy,  the  propagation  distance  was  about  0.6 
cm.  The  researchers  also  mention  that  energy  transport  by  means  of 
an  electron  beam,  while  taking  into  account  scattering  and  gas  heating, 
still  has  not  been  investigated  sufficiently  [57]. 


IX.  MFTI  RESEARCH  ON  BEAM  INSTABILITY 


Researchers  of  the  Moscow  Physico-Technical  Institute  (MFTI), 
under  the  leadership  of  K.  V.  Khodatayev,  carried  out  investigations 
from  1973  to  the  present  on  the  theory  of  the  resistive  hose  instability 
and  the  development  of  beam  oscillations  and  noise. 

A  nonlinear  solution  to  the  sausage  mode  (azimuthal  quantum 
number  m  =  0)  of  the  resistive  instability  was  determined  in  1984  [63] 
for  an  IREB  passing  through  a  conducting  plasma.  It  was  determined 
that  the  sausage  mode  of  the  instability  actually  does  develop.  In  the 
case  of  a  Bennett  distribution  function,  the  sausage  mode  of  the  resis¬ 
tive  instability  has  been  observed  to  be  dampened. 

Using  different  values  of  exitation  frequencies,  these  calculations 
showed  that  phase  mixing  leads  to  complete  damping  of  the  sausage 
mode  of  the  resistive  instability  in  the  case  of  a  hot  beam  with  a  Ben¬ 
nett  density  distribution.  According  to  MFTI  researchers,  this  conclu¬ 
sion  holds  true  in  conjunction  with  the  experimental  results  of  Refs.  64 
and  65  in  which  the  hose  instability  with  mode  m  =  1  was  observed, 
but  where  there  was  no  observation  of  the  sausage  mode  of  the  resis¬ 
tive  instability. 

MFTI  researchers  have  looked  at  the  energy  of  beam  noise  as  a 
function  of  time.  They  have  seen  this  beam  noise  mount  rapidly,  reach 
a  saturation,  and  then  slowly  decrease  with  time  [66].  They  have 
determined  that  the  beam  distribution  gradually  stops  varying,  the  high 
frequency  oscillations  become  dampened,  and  these  oscillations  transfer 
their  energy  to  the  thermal  movement  of  the  beam  particles.  Here  the 
beam  emittance,  which  is  the  characteristic  of  the  thermal  energy  of 
the  beam,  begins  to  increase. 

Equations  of  nonequilibrium  states  of  a  hot  ribbon  electron  beam 
moving  through  a  stationary  ion  background  indicate  that  the  beam 
relaxes  to  an  equilibrium  condition.  That  process  was  broken  down 
into  two  cases.  In  the  first  case,  the  transformation  took  place  in  con¬ 
junction  with  the  collisionless  theory  in  which  the  phase  mixing  of 
electron  trajectories  led  to  a  large  modulation  of  the  beam  distribution 
function.  The  calculations  showed  that  this  quenches  the  oscillations 
of  the  macroscopic  beam  characteristics.  In  the  second  case,  reversible 
processes  smooth  out  the  inhomogeneities  of  the  distribution  function 
in  phase  space,  and  thus  produce  beam  relaxation. 

The  characteristic  time  for  beam  relaxation  was  determined  to  be 
the  electron  oscillation  period.  The  rate  of  oscillation  damping  of  the 
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macroscopic  beam  characteristics  depends  upon  the  type  of  initial  exci¬ 
tation  of  the  beam  distribution.  Small  amplitude  oscillations,  which 
represent  a  fairly  small  area  of  the  phase  space,  are  damped  slower 
than  the  large  excitations,  which  possess  a  larger  phase  space  [66]. 


X.  IEM  RESEARCH 


Researchers  at  the  Institute  of  Experimental  Meteorology  (IEM)  in 
Obninsk  were  working  on  the  transport  of  IREBs  through  neutral 
gases  in  1980-1981  [67,68],  In  their  theoretical  studies,  they  have 
analyzed  a  condition  of  self-channeling  of  the  electron  beam  where  the 
action  of  the  beam’s  self-magnetic  field  is  capable  of  considerably 
decreasing  the  beam’s  spreading  due  to  multiple  scattering.  In  these 
investigations,  an  example  was  chosen  with  a  20  MeV,  34  kA  electron 
beam  transported  in  air  with  a  density  of  1.3  x  10”4  g/cm3  with  an  ini¬ 
tial  beam  cross  section  of  100  cm2.  It  was  determined  that  the  cross 
sectional  area  at  100  meters  increased  a  factor  of  4.5.  When  the  beam 
current  was  taken  at  low  values  ( I b  =  100  A),  however,  this  value 
increased  a  factor  of  3  for  all  other  conditions  unchanged  [67]. 

These  researchers  reached  the  following  conclusions  [67]: 

1.  The  condition  of  self-channeling  is  analogous  to  the  equilib¬ 
rium  case  for  the  pinch  effect  where  the  transverse  energy  of 
the  beam  particles  acts  as  the  temperature. 

2.  The  effect  of  the  beam’s  self-magnetic  field  can  considerably 
decrease  beam  spreading  caused  by  multiple  scattering. 

3.  Even  with  self-channeling,  the  beam  spreads;  initially  this  fol¬ 
lows  an  exponential  law. 

These  observations,  however,  have  been  well  documented  by 
Western  sources  in  the  early  1960s  and  1970s. 
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XI.  CONCLUSIONS 


Soviet  research  on  the  propagation  of  intense  relativistic  electron 
beams  appears  to  be  broad  and  systematic,  and  it  has  gathered  momen¬ 
tum  during  the  past  five  years.  One  obvious  application  of  such 
research  is  in  the  area  of  inertial  fusion  reactors  driven  by  electron 
beams.  Another  is  to  resolve  a  problem  in  the  quest  for  the  charged- 
particle  beam  weapon — the  feasibility  of  efficiently  transporting  such 
beams  through  atmospheric  air.  The  data  presented  here  and  drawn 
from  open  Soviet  publications  are  pertinent  to  some  basic  issues  of 
beam  weapon  research.  However,  they  cannot  be  directly  attributed  to 
such  research,  primarily  because  of  the  low  beam  energies  (<1.5  MeV), 
short  propagation  distances  (<  3  m),  and  low  pressures  used  in  the 
published  Soviet  experiments. 

An  indirect  hint  of  a  possible  military  involvement  is  the  disappear¬ 
ance  from  published  journals  of  five  key  experimental  teams  during  the 
period  1978  to  1981,  including  a  IYaF-TPI  team  lead  by  Yu.  P.  Usov 
and  two  teams  at  KhFTI  led  by  Ya.  B.  Faynberg,  Yu.  V.  Tkach,  and  A. 
K.  Berezin  (see  Fig.  1).  If  the  disappearance  of  these  teams  was  due  to 
security  classification,  then  the  total  Soviet  effort  in  this  research  area 
may  be  much  larger  than  appears  in  the  published  literature.  Even  in 
the  absence  of  publications  from  the  five  experimental  research  teams, 
the  overall  number  of  Soviet  published  papers  has  steeply  increased 
since  1983,  and  the  number  of  authors  has  shown  a  similar  increase 
since  1980  (see  Fig.  2). 

Five  new  Soviet  research  efforts,  initiated  in  the  late  1970s  or  early 
1980s,  are  currently  active  in  IREB  propagation  studies.  Of  these, 
FIAN-IOF,  led  by  A.  A.  Rukhadze,  demonstrates  the  largest  experi¬ 
mental  effort;  the  largest  theoretical  effort  is  being  carried  out  at 
IYaF-TPI  under  the  leadership  of  V.  P.  Grigor’yev. 

Because  the  research  being  carried  out  by  the  new  research  teams 
was  basically  the  same  as  that  of  the  five  old  teams,  two  hypotheses 
are  offered: 

1.  The  entire  research  program  has  been  transferred  to  the  new 
teams,  while  the  old  team  membership  has  been  disbanded. 

2.  The  new  teams  assumed  only  the  part  of  the  program  suitable 
for  open  publication,  and  the  old  teams  commenced  a  new, 
advanced  phase  of  research  subject  to  classification. 
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The  second  hypothesis  is  more  probable,  because  of  the  total  absence 
of  the  previous  authors  in  the  literature  of  the  1980s.  It  is  unlikely 
that  disbanded  team  members  would  altogether  cease  contributing  to 
scientific  research. 

The  specific  research  issues  pursued  in  the  published  literature  by 
the  former  and  the  new  Soviet  research  teams  on  IREB  propagation 
through  air  for  the  past  15  years  are  as  follows: 

•  The  formation  of  a  plasma  channel  created  by  an  IREB  propa¬ 
gating  through  background  air  and  the  investigation  of  the  effect 
of  beam  parameters  upon  the  plasma  channel  parameters  (and 
vice  versa). 

•  Determination  of  the  background  air  pressure  for  the  optimum 
transport  of  an  IREB.  Two  optimum  pressure  ranges  were  deter¬ 
mined: 

—  The  “ion  focused  regime” — 0.06  to  0.09  Torr. 

—  The  “low  pressure  window” — 1  Torr. 

•  The  current  enhancement  effect,  whereby  the  current  induced  in 
the  plasma  by  the  IREB  is  larger  than  the  original  beam  current. 

•  The  effect  of  the  resistive  hose  instability  upon  IREB  propaga¬ 
tion. 

These  key  issues  of  Soviet  research  match  the  main  interests  of  U.S. 
research  on  IREB  propagation  through  high  pressure  air. 

Soviet  researchers  found  that  an  IREB  can  be  transported  in  large 
volume,  large  diameter  insulating  channels,  and  through  free  space  by 
forming  two  channels  for  the  return  current  flow:  one  in  the  beam  and 
the  other  in  a  thin  cylindrical  layer  around  the  beam.  The  latter 
shorts  the  beam  current  to  the  accelerator  while  maintaining  beam 
self-focusing;  and  the  plasma  layer  has  a  stabilizing  effect  upon  the 
lateral  oscillations  of  the  beam,  thus  improving  IREB  transport.  This 
layer  was  assumed  to  be  a  plasma  channel  carrying  part  of  the  return 
current  of  the  beam.  When  the  beam  developed  hose  instability,  the 
cylindrical  symmetry  of  the  layer  was  not  disturbed,  but  the  entire 
beam  was  observed  to  be  displaced  laterally.  Measurements  revealed 
that  no  high  energy  electrons  moving  parallel  to  the  beam  flow  were 
present  in  this  layer  region.  It  was  determined  that  the  plasma  chan¬ 
nel  was  produced  by  the  ionization  of  the  background  gas  by  electrons 
of  the  IREB  leading  front,  and  the  beam  trajectory  was  determined  by 
the  movement  of  the  low  energy  electrons  and  not  by  the  main  part  of 
the  beam.  The  plasma  channel  remained  intact  longer  than  100  ns. 
The  beam  was  also  found  capable  of  tracking  the  plasma  channel. 

In  their  investigation  of  “current  enhancement,”  Soviet  researchers 
discovered  that  most  of  the  plasma  current  density  Jp  flowed  within 


64 


the  confines  of  the  electron  beam  in  a  time  about  100  ns,  which  is  con¬ 
siderably  longer  than  the  beam  pulse.  The  return  current  enhance¬ 
ment  (RCE)  effect  was  found  to  be  present  only  in  the  case  where  the 
IREB  had  a  large  energy  spread.  The  limitations  of  the  “two-stream 
instability”  at  the  low  pressure  of  IREB  propagation  has  been  reported 
only  by  NIIEFA,  and  the  earlier  Soviet  reports  have  mentioned  that  all 
current  enhancement  is  due  to  the  two-stream  instability.  However,  in 
1984,  FIAN-IOF  researchers  observed  a  current  enhancement  condition 
that  takes  place  at  pressures  close  to  atmospheric  when  none  of  the 
earlier  studied  causes  of  current  enhancement  could  be  the  reason  for 
the  effect.  This  anomalous  current  enhancement  was  found  to  increase 
with  pressure  and  propagation  length  and  was  in  the  direction  of  the 
IREB.  The  Soviet  researchers  claim  that  they  are  currently  making 
further  studies  of  this  effect. 

In  the  Soviet  open  literature,  the  only  mention  of  laser  guidance 
(the  use  of  a  laser  to  create  a  low-pressure  channel  along  whose  axis  an 
IREB  can  propagate)  was  made  in  1974  at  FIAN-IOF  by  G.  A. 
Askafyan  [69].  Here  a  preliminary  experiment  demonstrated  the  for¬ 
mation  by  a  laser  of  a  low-density  channel  in  air.  This  channel  had  a 
finite  lifetime  of  several  hundreds  of  microseconds  with  a  density  10  to 
20  times  lower  inside  the  channel.  A  50  keV  electron  beam  pulse  was 
passed  through  this  channel  successfully,  thus  demonstrating  laser 
guidance.  Other  FIAN-IOF  papers  published  by  G.  A.  Askar’yan  and 
his  co-workers  have  made  statements  on  the  use  of  laser  guidance  for 
the  successful  transport  of  charged  particle  beams.  There  are  no  recent 
publications  on  laser  guidance,  and  it  is  believed  that  this  subject  is 
censored  by  the  Soviet  authorities. 

Most  of  the  Soviet  reports  that  have  been  considered  here  do  not 
specifically  mention  a  single  application  of  their  research.  However, 
many  papers  generally  indicate  the  importance  of  optimum  propagation 
conditions  for  high  energy,  high  current  beams,  and  their  transport  to  a 
maximum  allowable  distance.  In  a  few  cases  (VEI,  NIITP,  MFTI), 
where  the  injection  of  the  beam  into  atmospheric  pressure  air  was  con¬ 
sidered,  specific  applications  for  beam  welding,  beam  cutting,  and 
specimen  hole  boring  were  mentioned  along  with  the  general  objective 
of  achieving  maximum  propagation  length  in  air. 


Appendix 


SOVIET  RESEARCH  TEAMS  INVOLVED  IN 
IREB  PROPAGATION 


Soviet  research  teams  involved  in  IREB  propagation  in  air  at  pres¬ 
sures  from  0.1  to  760  Torr  are  listed  below  by  institute  affiliation. 
Group  leaders  are  indicated  in  italics. 

IYaF-TPI:  Institute  of  Nuclear  Physics  at  Tomsk  Polytechnic  Institute,  Tomsk 
Team  1:  Theory  and  Computer  Simulation 
V.  P.  Grigor’yev 

A.  G.  Potashev  N.  S.  Shulayev  A.  V.  Zakharov 

A.  I.  Koryakin  A.  S.  Popov  A.  N.  Didenko 

G.  P.  Isayev 

Team  2:  Experimental 
Yu.  P.  Usov 

A.  N.  Didenko  V.  A.  Tuzov  A.  I.  Ryabchikov 

Ya.  Ye,  Krasik  A.  V.  Petrov  G.  I.  Kotlyarevskiy 

IAE:  Kurchatov  Institute  of  Atomic  Energy,  Moscow 
L.  I.  Rudakov 

V.  P.  Smirnov  V.  V.  Gorev  S.  V.  Zakharov 

L.  Ye.  ATanchuk  V.  D.  Vikharev  S.  F.  Zhandarov 

V.  D.  Korolev  L.  I.  Urutskoyev  S.  F.  Grigor’yev 

A.  V.  Gordeyev  V.  V.  Zazhivikhin  V.  Yu.  Shuvayev 

FLAN:  Lebedev  Physics  Institute/IOF:  Institute  of  General  Physics,  Moscow 
A.  A.  Rukkadze 

A.  A.  Savin  G.  P.  Mkheidze  A.  L.  Ipatov 

G.  A.  Askar’yan  B.  M.  Manzon  R.  R.  Kikvidze 

I.  M.  Minayev  A.  G.  Shkvarunets  V.  Yu.  Shafer 

S.  N.  Kabanov  A.  A.  Korolev  A.  I.  Kuz’min 

V.  Ye.  KuTbeda  V.  Yu.  Nosach 

KhFTI:  Physical-Technical  Institute,  Khar’kov 
Team  1 

Ya.  B.  Faynberg 
Yu.  V.  Tkach 

1. 1.  Magda  G.  V.  Skachek  S.  S.  Pushkarev 
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V.  A.  Bondarenko 
S.  I.  Naysteter 
P.  M.  Lebedev 

Team  2 

Ya.  B.  Faynberg 

A.  K.  Berezin 
V.  A.  Kiselev 

B.  G.  Safronov 

N.  M.  Zemlyanskiy 

L.  I.  Bolotin 
V.  A.  Buts 


I.  P.  Panchenko 
I.  N.  Onishchenko 


V.  P.  Zeydlits 
G.  L.  Fursov 
I.  K.  NikoFskiy 
S,  S.  Krivulya 
V.  I.  Kurilko 


B.  D.  Panasenko 
V.  I.  Shevchenko 


I.  A.  Grishayev 
I.  N.  Onishchenko 
V.  I.  Mirnyy 
A.  M.  Yegorov 
A.  P.  Tolstoluzhskiy 


NIIEFA:  Yefremov  Institute  of  Electrophysical  Equipment,  Leningrad 

L.  V.  Dubovoy 

T.  S.  Gosteva  G.  R.  Zablotskaya  B.  A.  Ivanov 

V.  I.  Chernobrovin  S.  A.  Kolyubakin 

NIITP:  Scientific  Research  Center  of  Thermal  Processes,  Moscow 
V.  M.  Iyeulev 

A.  S.  Koroteyev  V.  V.  Koba  I.  G.  Kulakov 

VEI:  All-Union  Electric  Engineering  Institute,  Moscow 

A.  V.  Zharinov 

M.  A.  Vlasov 

S.  V.  Nikonov  S.  I.  Vybornov  V.  A.  Safonov 

A.  V.  Rykhlov  V.  P.  Popovich  Yu.  A.  Yel’chin 

V.  A.  Artemov  L.  G.  Dubas  I.  P.  Denisova 

A.  V.  Smirnov 


MFTI:  Moscow  Physico -Technical  Institute,  Moscow 
K.  V.  Khodatayev 

V.  G.  Leyman  A.  A.  Ivanov  K.  G.  Gureyev 

V.  0.  Zolotarev  S.  D.  Stolbetsov  N.  Ye.  Rozanov 

N.  S.  Putvinskaya  0.  B.  Ovsyanikova  F.  F.  Kamenets 

IEM:  Institute  of  Experimental  Meteorology,  Obninsk 

A.  P .  Budnik 
P.  N.  Svirkunov 
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